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5.^S®[EGIE - MELLON UNIVE. RSITY 
1969-1970 
PART I. GENERAL REPORT 


A . Introduction 


This report describes research in the areas of nuclear chemistry, 
geochemistry, and cosmochemistry carried out on the main campus of Carnegie - 
Mellon University from summer 19&9 through spring 1970 under the principal 
sponsorship of the United States Atomic Energy Commission under Contract No. 
AT(30-l)-844. 

We include in this document also reports on the progress of research 
sponsored by the National Aeronautics and Space Administration under its Lunar 
Samples Analysis Program under Contract No. NAS -9 -8073 for "Radiogenic Isotope 
Studies of Lunar Materials" . From the viewpoint of nuclear chemistry, this work 
is an extension of studies of nuclear phenomena in nature utilizing terrestrial 
and particularly meteoritic materials, and an opportunity to apply analytical 
and interpretative techniques to some extent already developed to a new and 
important class of materials . 

The principal part of this Report is a series of contributions authored 
by the students, research associates, and other research personnel who have 
participated in one phase of the work or another. Some of the contributions 
describe incomplete work, plans, ideas, and calculations. Many of the numerical 
data and suggested interpretations presented here are preliminary or tentative. 
Those wishing to use or quote such information are requested to treat it as such, 
and if possible to check as to whether final results have been published or 
prepared for publication. We will be happy to supply additional information 
through correspondence, preprints, or reprints as available. 

B. Personnel 


1. Principal Investigator 

Truman P. Kohrnan: A.B. (1938) Harvard College; Ph.D. (1943), University 

of Wisconsin; Professor of Chemistry. 

2. Research Associates 

Dr. Lance P. Black: B.S. (1966) The Australian National University; 

Ph.D. (1969) The Australian National University; Research Geochemist and Co-Investigator, 
Lunar Sample Analysis Program. 

Dr Haruhiko Ihochi: B.S. (1963) St. Paul's University, Tokyo; Ph.D. ( 1969 ) 

University of Arkansas; Research Chemist and Co-Investigator, Lunar Sample Analysis 
Program. 



Dr. Brio M. P. Trivedi: B.Sc. (1962) Indian Institute of Technology, Kanpur ^ 

India j Ph.D. (1969) Indian Institute of Technology, Kanpur, India; Research Chemists 

3 . Graduate Students 


Mr. James M. Huey: B.S. (1966) University of Kentucky? M.S. (1969) 

Carnegie Mellon University; fourth-year graduate student in Chemistry. For the past 
three years Mr* Huey has been a NASA. Graduate Trainee. 

Mr. Ramzi Y. Saleh: B.S. (1963) Ain Shams University, Cairo, Egypt, U.A.R.J 

second -year graduate student in Chemistry until May 1970* (See below). 

Mrs. Nancy Man-Na (nee Kuo) Kan: B.S. (1969) National Taiwan University? 

first -year graduate student in Chemistry. 

4. Undergraduate Students 


Miss Nina A. Landes: Chemistry, ' 70? full-time during summer 1969, 

part-time during academic year. (See below). 

Mr. Charles J. Vukotich: Chemistry -Psychology . * 71? full-time during summer 

1969, part-time during academic year. 

Miss Sarah E. Weber: Chemistry, ! 71? part-time during academic year. 

Mr. Mark J. Yeager: Chemistry, *71; full-time during summer 1 969, part-time 

during academic year and summer 1970* ' 

Mr. John R. Yocum: Chemistry, 1 73 ? part-time during academic year and summer 

1970. . 

3. Technical Assistants 

Mr. Mark W. Haramic: Diploma (1968) Air Force Electronics Course for 

Aircraft Radio; Diploma (l 969) Penn Technical Institute; Research Electronics 
Technician. 

.Miss Nina A. Landes: B.S. in Chemistry (1970) Carnegie -Mellon University; 

Junior Research Chemist, summer 1970 (see above). 

Mr. Ramzi Y. Saleh: (see above); M.S. in Chemistry (1970) Carnegie -Mellon 

University? Junior Research Chemist, summer 1970* 

Mr. J. David Ulmer: B.S. in Physics (1970) Carnegie -Mellon University; 

Junior Research Physicist, summer 1970. 

Mrs. Ignez de Carvalho: B.S. (196l) Universidade Federal de Minas Gerais, 

Minas Gerais, Brazil; Instructor, Universidade Federal de Minas Gerais (1964-1968)? 
Assistant Professor of Chemistry, Universidade Federal de Minas Gerais (1968); 
Laboratory Technician. 
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6. Won -Technical Assistants 


Mrs. Gail R. Chasey: Secretary through April 1970* 

Mrs . Judith W. Watt: Secretary beginning May 1970. 

C . Facilities 

1, Laboratories 

All of the work described herein has been conducted in the laboratories 
in Doherty Hall on the main campus which have been described in previous reports. 

These include three chemical laboratories and an electronics laboratory on the 
third floor; a radiochemical laboratory, parts of several chemical laboratories, 
and a clean room on the second floor; and a large air-conditioned room in the 
basement which houses the mass spectrometer and the low-level counting instruments. 
Although this space has served us well, it is now remote from the main focus of 
advanced study and research of the Department of Chemistry, which has now shifted 
to the Mellon Institute Building. Moreover, it has become increasingly apparent 
that the clean room, with ordinary air-conditioning and air filters, is inadequate 
for work with extremely low levels of trace elements, particularly lead. 

The space allocated to this research group in the Mellon Institute Building 
is now undergoing renovauion and preparation for all aspects of our work. It 
consists of the following units, all but the last of which are on the northeast 
corner of the fifth floor of the building: 

An office suite for professor, research associates, and secretary, including 
a library and reading room which can be used also for conferences. 

An electronics laboratory with a small instruments room. 

A large conventional chemical laboratory. 

A radiochemical suite consisting of a “hot" and a "warm" area. 

A clean laboratory, within which are a walk-in hood and a fume hood each 
with HEPA (High -Efficiency Particulate -Air) filters. 

On the north side of the first floor (basement), a large nuclear measurements 
laboratory, which will house the mass spectrometer and most of the counting instruments 
and nuclear spectrometers. 

All of this space is air-conditioned, the clean laboratory being provided with 
specially -filtered room air in addition to the HEPA filters. According to present 
plans, the renovations will be completed during July and the move will be made 
during August. 

2 . Instruments 

The radiation detection instruments available are essentially the same as 
described in last year's Report. Their present status is described in Section II. 

D. 1. of this Report. 
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The most important improvement during the year has been in the signal -handling 
part of the mass spectrometer, which has been equipped with a digital voltmeter, 
digital printer, and rapid beam-switching devicej see Section IX. D. 5» of this 
Report . 

D. Publications 

I i 

1, Published and in Press 

ISOTOPIC COMPOSITION OF METEORITIC THALLIUM 

R. G. Ostic, H. M. El-Badry, and T. P. Kohman 

Earth and Planetary Science Letters 7, 72-7 6 (1969 Oct.) 

NYC -844 -78 

LEAL AND THALLIUM ISOTOPES IN MARE TRANQUILLITATIS SURFACE MATERIAL 
T. P. Kohman, L. P. Black, H. Ihochi, and J. M. Huey 

Science 167 , 481-485 (1970 Jan. 50 )', Geochimica et Cosmochimica Acta 54, 
in press (1970 Supplement) 

CMU-MSA-21-5 (Revised) 

LEAST -SQUARES FITTING OF LATA WITH LARGE ERRORS 
T. P. Kohman 

Journal of Chemical Education 47, in press (197° July ?) 

NYO-844-77 

TABLES OF URANIUM -THORIUM -LEAD DECAY -GROWTH FUNCTIONS 
T. P. Kohman ' 

NYO-844-79 

2. Presented at National Scientific Meetings 

LEAD AND THALLIUM ISOTOPES IN MAKE TRANQUILLITATIS SURFACE MATERIAL 

T. P. Kohman, L. P. Black, H. Ihochi, and J. M. Huey 

Apollo 11 Lunar Science Conference, Houston, Texas, 1970 January 8 . 

CMU -NASA -21 -5 

A UNIVERSAL MAGNITUDE SCALE FOR ASTRONOMICAL OBJECTS 
Truman P. Kohman 

American Astronomical Society Meeting, Boulder, Colorado, 1970 June 10. 

NY0-844-80 

5 . In Preparation 

EXTINCT NATURAL RADIOACTIVITY AND POSSIBLE LUNAR MANIFESTATIONS 
T. P. Kohman 

A UNIVERSAL MAGNITUDE SYSTEM FOR ASTRONOMICAL OBJECTS 
T. P. Kohman 

CHRONOLOGY OF GALACTIC HEAVY-ELEMENT NUCLEOSYNTHESIS 
J. M. Huey and T, P. Kohman 

RADIONUCLIDE DISTRIBUTIONS IN THICK IRON TARGETS WITH 0.1 TO 6.1-GeV PROTONS 
G.V.S. Rayudu and T.P. Kohman 
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RADIONUCLIDE DISTRIBUTIONS IN THICK SILICATE TARGETS BOMBARDED BY 0.1-, OA- 
1.0-, AND 3.0-GeV PROTONS 
J. P. Shedlovsky and G. V. S. Rayudu 


u 


fc* 


rv. r* g~ xf* i/f 

MERCURY AND BISMUTH ^ABUIMCE^'IN rf«3N*METEORITES BY NEUTRON ACTIVATION 
J. T. Tanner and T. P. Kohman 


MERCURY, THALLIUM, AND BISMUTH IN THE METAL AND TROILITE PHASES OF IRON 
METEORITES 

S. N. Tandon and T. P. Kohman 

COSMOGENIC RADIONUCLIDES IN STONE METEORITES 
P. J. Cressy, Jr. and T. P. Kohman 

COSMOGENIC X-RAY AND p-RAY EMITTERS IN IRON METEORITES 
J. H. Kaye, M. M. Chakrabartty, and T. P. Kohman 

THORIUM ISOTOPES METHOD FOR DATING MARINE SEDIMENTS 
I. Almodovar and T. P. Kohman 

DATING OF MARINE SEDIMENTS BY IONIUM AND PROTACTINIUM METHODS 

T. P. Sarma and T. P. Kohman 
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during the past year, whether specifically mentioned in the individual sections of 
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professional colleagues who have sent us reports, preprints, and reprints of their 
publications . 



6 . 


PART II . RESEARCH PROGRESS REPORT AND PLANS 


A. Nuclear Geochemistry 


N71-10305 

II. A. 1 ■ - 


INVESTIGATION OF Sm 146 - 


Nd 142 AND Sm 147 


- Nd 143 AS DATING METHODS 


Ramzi Y. Saleh 


I. Introduction 


Investigation of Sm 1 ^-Ncf 1 ’ 42 and Sm 1 ^-Nd -1 43 . systems for 
their potential applications to geochronology (1,2,3) has been continued. 

The utilization of these systems as dating methods will depend 
on the extent of Sm-Nd fractionations ' in natural materials . Sm-Nd frac- 
tionations in various types of natural materials have been discussed 
previously (3) . 

The discussion of the ion-exchange separation of neodymium and 
samarium from each other as well as from the other rare-earth elements 
has been reported previously (3) . 

Mass spectrometry of 1-5-pg samples of standard solutions of 
natural neodymium and samarium (4) has been found to be quite feasible. 


II. Chemical Procedure 

Dissolution of rock samples and group separation of the rare- 
earth elements has been described elsewhere (3) . 

The separation of samarium and neodymium from each other as well 
as from other elements using 1-M ammonium lactate as eluant, Dowex 50W- 
X8 50-100 mesh as a resin, and 67°C as an operating temperature (1,2, 3, 6) 
has been found to be difficult due to the fact that hot columns require 
specialized 'equipment and techniques. Moreover, the separability of 
these two cations, which is a function of the ratio of their distribution 
coefficients, was not of the right order of magnitude for the' elution (7) . 

Another technique was used in this laboratory for the separation 
of samarium and neodymium (5) . Bio Rad AG-50W-X4 200-400 mesh was used 
as a resin. 0.4-M ammonium a-hydroxy isobutyrate was the eluant used. 

The design of the column was the same as that used for the 1-M ammonium 
lactate separation (3) but different dimensions (30-cm height and 0.6-cm 
diameter) . A coarse fritted-glass disc was sealed at the bottom of the 
column to support this small-particle size resin instead of using glass 
wool. The cation-exchange elution curve of samarium and neodymium was 
determined using two different methods. 

A, Determination of Sm and Nd Peak Positions by Visual Inspection 
of Their Oxalate Precipitates and Pm i4/ Tracer Peak by Beta- 
Counting 

The rare-earth preparation consisted of a mixture of 5 mg 
Sm, 5 mg Nd„ and Pm 147 tracer. A solution of 0.4-M ammonium a-hydroxy 
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isobutyrate of pH 4.2 was used 'as eluant. Ten ml was the. fraction size col- 
lected at room temperature. The resulting .elution curve showed that 
there was no separation at all, which can he attributed to the high pH 
value (4.2). 

The same experiment has been repeated, using the -same concentration 
of the rare-earth elements and the same technique, with the exception that 
the pH was reduced to 3.8. The resulting data showed incomplete separation, 
which might be. explained by the same reason given above. 

The separation of Sm, Pm 147 , and Nd was improved when the pH was 
reduced to a lower value of 3.6. Figure 1 is an ion-exchange elution 
curve of macroscopic amounts of Sm, Nd, and Pm 147 tracer. The resulting 
elution curve obtained using the technique given above showed that a 
sharp separation between the three rare-earth peaks has been obtained. 

The results of this experiment were confirmed when the same experiment, 
with the same technique and conditions, was repeated-. 

B. Tracer Experiments 

, The next preparation consisted of a mixture of Sm I4S , Pm 145 , and 
Nd 147 gamma-emitting tracers. One mg La was added as a carrier. This 
experiment was done at room temperature, and the flow rate was 3 drops/ 
min. The same technique was repeated with the exception that the volume 
of the fraction collected was 5 ml each. Gamma counting was done on 
every 5-ml test tube, using a Ge(Li) gamma-ray detector and a. 512-channel 
analyzer. Sm 145 was detected by 61-keV y-ray.s, Pm 145 by. 72-keV x-rays, 
and Nd 147 by 91-keV y-rays. Figure 2 is a semi-logarithmic plot of 
counting rate versus the volume of the eluate. This figure shows that a 
sharp separation has been obtained. 

Another tracer experiment has been done using a rare-earth mineral. 

A 0.5-g sample of allanite was decomposed using the chemical procedure 
described elsewhere (1,3). The total rare-earth elements were separated 
as a mixture, dissolved .in HC1, and kept as a stock solution. A mixture 
of one- tenth of this stock solution, which is equivalent to 50-mg allanite, 
and which contained about 10-mg rare-earth elements, and tracers of Sm 145 , 
Pm 145 , and Nd 147 was prepared. The previous conditions for. the Sm and Nd 
separation were .repeated. The data obtained are shown in Figure 3-a,. It 
is observed from this figure that Sm, Pm, and Nd were eluted rather more 
rapidly than the corresponding peaks in Figure 2. This can be attributed 
to the fact that in the present experiment there was an appreciable fraction 
of the column -volume nearly saturated by the rare-earth content-. 

The rare-ear-th concentration, was reduced, when a mixture of 5-mg ' 
allanite (one hundredth of the stock solution was used which contained 
1000 Ug rare-earth elements) and tracers of Sm 145 , Pm 145 , and Nd 147 
(containing ~ 15 pg rare-earth elements) was loaded on the column. Figure 
3-b'is the resulting elution curve. 

The concentration of the rare-earth in. the ion-exchange column 
was still further reduced so as to determine the concentration range in 
which the Sm, Pm, and Nd .peak positions will be independent of the amount 
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present. A mixture of 0.5-mg allanite (containing approximately 100 yg 
rare-earth elements) and the three radioactive tracers (containing ~ 15 yg 
rare-earth elements) were loaded on the column. The resulting elution 
curve is shown in Figure 3-c. Moreover, one more experiment was done 
using Sm 145 , Pm 145 , and Nd 147 tracers. The total rare-earth element 
content in the column was ~ 15 yg. The data of this experiment which are 
represented in Figure 3-d were found to be close to those in Figure 3-c. 
These close results obtained (Figures 3-c and 3-d) may result in the 
conclusion that if the rare-earth concentration in the ion-exchange column 
is 0.1 mg, the Sm, Pm, and Nd peaks will be located at the same positions 
in the elution curve. 

III. Mass Spectroscopy 

The technique of mass spectrometry of samarium and neodymium 
is described elsewhere in this report (4) . 


IV. Isotope Dilution Analysis of Sm and Nd 

Isotope dilution technique will be used to determine the concen- 
trations of samarium and neodymium in natural materials. 

Samarium spike enriched in 150 (Sm 150 fraction and Sm 1 5 0 / Sm 1 4 9 ratio 
as provided from the analysis sheet were 95.1% and 80 .59 > respectively) 
and neodymium enriched in 150 (Nd 150 atom and Nd 150 /Nd llf6 ratio were 94.7 5% 
and 86. 14, respectively) were obtained from Oak Ridge National Laboratory, 

Oak Ridge, Tennessee. Each spike was dissolved in 10 ml d'il. nitric acid 
kept in a clean 10-ml volumetric flask, and designated as Sm 150 spike, 

1970 and Nd 15 ' 0 spike,. 1970. 

Standard nitrate solutions of reagent samarium and neodymium 
oxides of accurately known concentrations were prepared for use in cali- 
bration of the spikes. A very sensitive Sartorius Balance was used to 
determine the accurate weights of these reagents. These nitrate solutions 
will be standardized using gravimetric techniques, since the stoichiometric 
composition has not been verified. These reagents were manufactured by 
Alfa Inorganics, Inc., Beverly, Massachusetts. The concentrations of 
these nitrate solutions were 77 yg Sm/10 yl and 39.99 yg Nd/10 yl and 
labeled as Sm-STD. S0LN. 1970 and Nd-STD. SOLN. 1970. 

A mixture of 10 yl Sm 150 , spike solution (~ 4.3 yg Sm 150 ) and 10 yl 
of the standard reagent (77.00 yg natural Sm) estimated to yield a Sm 150 / 
Sm 149 ratio of unity was prepared and analyzed by mass spectrometry (4). 
Another mixture of 10 yl Nd 1 0 spike solution (~ 4.1 yg Nd 1 - 30 ) and a 
10 yl of natural neodymium solution, (39.99 yg natural Nd) was prepared 
and analyzed in a similar manner. 

Table I shows the different measured isotope ratio of samarium 
and neodymium and their mean standard deviations ( 0 ) . Sm 149 and Nd 146 

were selected as reference isotopes for the calculations of the Sm and 
Nd ratios, respectively. The exact concentrations of both Sm 150 and Nd 150 
spike solutions were 4.35 yg Sm is °/10 yl and 4.11 yg Nd lso /10 yl, respec- 
tively. The Sm 1 50 /Sm 1 49 and Nd 150 /Nd 146 ratios used for this calculation 
were 0.5340 and 0.3291, respectively (4). 
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About 1 g of very old materials (Precambrian rocks and minerals, 
meteorites, and lunar samples) will be chemically decomposed according 
to the procedure described elsewhere (3) . This solution will be 
divided into two fractions. One of these fractions will be spiked for 
the determination of Sm and Nd concentrations. The other fraction will 
be used for the determination of the- isotopic composition of neodymium. 
It is conceivable that very careful neodymium measurements in very old 
materials might reveal Nd 143 and Nd 142 abundance variations which may 
be useful for chronological interpretations, particularly if attention 
is paid to variation in the, Sm/Nd ratio. Each of these two fractions 
will be subjected to the procedure described in part II. B of this 
section. 


References 

(1) R. Y. Saleh, Semi-Annual Progress Report to NASA, Contract NAS-9-8073 

(1 May, 1968 - 31 October, 1968), Report CMU-NASA-21-1 , Section VI. 

(2) R. Y. Saleh, Semi-Annual Progress Report to NASA, Contract NAS-9-8073 

(15 May, 1969), Report CMU-NASA-21-2 , Section VI. 

(3) R. Y. Saleh, 1968-1969 Progress Report, USAEC Report NYO-844-76 , 

II. B. 2. 

(4) L. P. Black, This Report, II. D. 3. 

(5) H. L. Smith and D. C. Hoffman, J. Inorg. Chem. , 3_, 243-247 (1956). 

(6) R. A. Schmitt and R. H. Smith, "A Program of Research for the 

Determination of Rare-Earth Elemental Abundances in Meteorites", 
General Atomic Report 3411 , September 15, 1961 through August 
14, 1962. 

(7) 0. Samuelson, Ion Exchange Separation in Analytical Chemistry , 

(John Wiley & Sons, Inc., New York, 1963) Chapter 3. 
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TABLE I 

MASS- SPECTROMETRIC BATA FOR SPIKE CALIBRATION 


* ■ - 

Samarium 

Neodymium 

Ratio 

Mean Value 

a of Mean 

Ratio 

Mean Value 

a 

of Mean 

Sm 144 /Sm 149 

0.2277 

± 

0.0002 

Nd 142 /Nd 146 

1.5941 

± 

0.0076 

Sm 147 /Sm 149 

1.0884 

± 

0.0008 

Nd 143 /Nd 146 

' 0.7104 

+ 

0.0004 

Sm 14s /Sm 149 

0.8157 

+ 

0.0005 

Nd I44 /Nd 146 

1.3879 

+ 

0.0008 

Sm 1 50 /Sm 1 49 

1.0841 

+ 

0.0008 

Nd 145 /Nd 145 

0.4837 

+ 

0.0009 

Sm 1 52 /Sm 1 49 

0.5212 

± 

0.0005 

Nd 148 /Nd 146 

2.9665 

+ 

0.0087 

Sm 154 /Sm 149 

0.6197 ’ 

i 

+ 

0.0004 

Nd 150 /Nd 146 

1.1545 

± 

0.0003 
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Figure 1. Cation-exchange separation of macroscopic amounts of Sm and Nd 
and tracer amounts of Pm 147 . (Exchanger, AG 50W-X4, 200-400 
mesh; column, 30 cm long by 0.6 cm diameter; eluant, 0. 4-M 
ammonium a-hydroxy isobutyrate, pH = 3.6; flow rate = 3 drops/ 
min; T = room temperature; each test tube contains 10 ml). 

(a) 5 mg Sm. (b) Pm 147 tracer, (c) 5 mg Nd. 


Relative Size of Sm and Nd Oxalate Precipitates 



Counting Rate (cpm) 
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TABLES OF UHANIUM -THORIUM -LEAD DECAY -GROWTH FUNCTIONS 

Truman P. Kohman 


Introduction 


Calculations involving the tf^-Pb 206 , U 235 ^ 207 , and Th 232 -Pb 20S 
radioactive systems, particularly in geochronology applications, involve 
the use of exponential decay or growth functions, singly or in combinations , ' ’ 
of two or four exponentials involving a single time or two different times. ' 
Often the times are the quantities to be evaluated by means of these functions, 
but the combination functions cannot be solved analytically for the times. 
Consequently the solutions are usually obtained by successive trials or 
graphical methods, requiring that the functions be computed several times. 

To eliminate some of the tedium of these calculations, a set of six 
tables giving a total of approximately 15,000 values of functions have been 
computed and printed with an electronic computer. 

Functions Computed 


The following quantities are parameters or arguments of the various functions 
which were computed: 

A g = disintegration constant of U 233 

Xs - disintegration constant of U 238 

X 2 = disintegration constant of Th 232 

C = (U 238 /^ 35 ) (atomic ratio) 

present} 

« - l/c - (u 233 /u 238 )p resent (atomic ratio) 

T = age of beginning of single-stage interval ending at present 

Ti = age of end of multiple -stage interval in past 

T a = age of beginning of multiple -stage interval in past 

The single -stage -model functions tabulated in Table I are: 


X S T 

e -A 5 T 

a 8 t 

e "AgT 

x 2 t 

-A 2 T 

e 
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rf A r— - rp ^ 

. a 5 t 

A 5 e 5 

JT( \ 

6 

a 2 t 

a 5j ) 

" “ x e A a T 

a 8 e 

A 2> A Q* 1 / 

As e 

A S T 


a A = T - 1 


e XsT 

- 1 

s(«,a 5 , a 8 ,t) 


v(a 2 ,a 8 ,t) 

= T5T 

- 1 


The multiple -stage functions tabulated are: 


Table IX: 

e A 5 T 2 _ e A sTl 



Table III: 

e A 8 T 2 _ e ^sTi 



Table IV:. 

e A 2 T 2 _ e ^2^i 

- 



* 

e AsT 2 ■ 

. e ^ Tl 

Table V: 

S( (X, A 8 , A 8 ,T 2 ,Ti ) 

" “ 

- e AsTl 



e A 2 T 2 

- e AaTl 

Table VI: 

V(.A 2 ,A 8 ,T 2 ,T3.) 

e AaT2 

. e A8 ^ 


The physical significance of some of these functions is as follows: 

R(cc,A 5 ,A 8 ,T) - instantaneous production -rate ratio of Pb 207 to 

Pb 206 in natural uranium at time -T (present time = 0). 

S(a,A 5 ,A 8 ,T) . = slope of single-stage isochron for interval from -T to 

present on Pb 207 /Pb 204 versus Pb 206 /Pb a ° 4 diagram. 

U( A 2 ,A a ,T) = instantaneous production-rate ratio of Pb 208 to Fb^ cs 

at time -T in system with equal numbers of Th 232 and 
U 238 atoms at present. 

V(A 2 ,A 8 ,T) ■ = factor which, when multiplied by present Th 23a /U 23S 

ratio, gives slope of single-stage isochron for 
interval from -T to present on Pb 2 ° 8 /Pb 204 versus 
Pb 206 /Pb 204 diagram. 

S(ex, a 5 ,A 3 ,T 2 ,Tx) — slope of multiple -stage isochron for interval from -T 2 

to -a?! on Ph 2 ° 7 /Pb 204 versus Pb 206 /?! 204 diagram. 

= factor which, when multiplied by present Th 232 /U 238 ratio, 
gives slope of multiple -stage isochron for interval from 
. -T 2 to -Ti on Pb 208 /Pb 204 versus Pb 20S /Ph 204 diagram. 


V(A 2 ,As,T 2 ,Ti ) 
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Nuclear Data Used 


The computations have been made for the values of the nuclear constants 
adopted by Russell and Farquhar (l): 

\ 5 = 0.9722 x 10 “ 9 y -1 

A s = ■ 0.1537 X 10 " 9 y' 1 

A a = 0.0499 X 10 " 9 y _1 

C = 137.8 

These are very close to the earliest reasonably accurate values (2, 3 .>4, 5*6). 

There is some evidence that more accurate values are available (1,7,8). However, 
these are probably not far from the true values, and they have been recommended 
again recently by Kanasewich (9). Since most workers in lead isotope geochronology 
use these values, it seems best to retain them for all calculations until agreement 
can be reached on a better set. 

Description of Tables 

The single-stage functions are evaluated at intervals of 0.01 Gy = 10 My 
for ages up to 6 Gy = 6,000 My. Linear interpolation should be quite adequate 
for age intervals of 0.001 Gy = 1 My. 

The multiple -stage functions are evaluated for T 2 > Ti at intervals of each 
argument of 0.1 Gy = 100 My for the same range. For 0.01 -Gy (10 -My) intervals 
linear interpolation should suffice except for functions involving A5 and Ti ~ T 2 . 
In such cases, and in all cases where accurate values are required for 0.001 -Gy 
(l-My) intervals. Table I can be used with an intermediate amount of tedium. 
Alternatively, tabulated values in the region of interest can be plotted and 
connected by smooth curves, and intermediate values can he determined by graphical 
interpolation. 

The computations were made by the Carnegie -Mellon University Computation 
Center’s UNIVAC 1108 computer with a FORTRAN V program designated URTHLE. The 
exponentials in the first three columns of Table I were evaluated directly, and 
all other values in Tables I -VI were computed from them. Appropriate derivative 
functions were substituted for limiting cases which would lead to o/o. Compilation 
and execution required 23 seconds . 

Availability of Tables 


These tables are available as CMU-AEC Report No. NYO -844-79, dated 1970 May 15 . 
A listing of the program can be supplied on request. 
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II. A. 3. 

SEARCH FOR LONG-LIVED TRANSURANIUM ALPHA EMITTERS AND OTHER ANOMALOUS 
ALPHA EMITTERS IN NATURAL MATERIALS 

Nina Landes 


I. Introduction 


Any discussion of a search for unknown transuranium elements in 
natural materials must necessarily begin with a theoretical justification 
for their survival as possible primary alpha emitters. First it must be 
shown that superheavy elements could have been formed in nucleosynthesis. 

Then, it must be proven that the elements thus formed could not only have 
outlasted the period At when the elements reacted, fractionated, and con- 
densed into the minerals composing the planets and meteorites, thus quali- 
fying them to be classified as "extinct natural radioactivity" (1) , but 
also be around today. Nilsson et al. , point out that t 1/s must be greater 
than 2 x 10 s years to be detected by the techniques previously tried by 
Thompson, by Ghiorso, and by Fowler (2). 

Superheavy elements could have been produced in the r- process during 
a Type I supernova explosion (3), as evidently were considerable portions 
of the known heavy elements, besides all the isotopes of uranium and thorium, 
and Pu 244 . Difficulties arise because the heavier the intermediates formed 
by neutron bombardment, the more likely they are to fission spontaneously 
before they can beta decay to stable superheavy nuclides. Still, the pos- 
sibility that some would have survived cannot be ruled out. Once formed, 
the nuclide must be sufficiently stable to alpha and beta decay and to spon- 
taneous fission, the latter tending to assume more and the second less 
importance with increasing mass number in the transuranium region (3). An 
account of the projected half lives for all three modes led Nilsson (2) to 
propose that the nuclide 294 110 is likely to be around today. Large uncer- 
tainties in the calculations give leeway for other elements instead of 294 110. 
Seaborg concluded that it is feasible to expect the survival of a small number 
of primary emitters, particularly in the islands of stability around Z = 126, 

N = 184, and Z = 114, N = 184 (3). 

Studies have been made of the newly discovered phenomenon of shape 
isomerism, in which a relative minimum can exist in' the potentials of known 
elements due to a stable vibrational state conductive to fission (4). This 
opens up the possibility of various modes of decay from this state in the 
known elements, which have previously been predicted to be energetically un- 
likely, and poses a rationale for anomalous alpha activity. 

Many investigators have reported evidences of what they interpret 
as extinct natural' radioactivity, some of which might otherwise be inter- 
preted as being caused by previously unknown primary radioactive emitters. 
Observances of an unusual abundance of Xe 129 have been interpreted as evi- 
dence of. "in situ" beta decay of an extinct I 129 parent, which is itself a 
daughter of spontaneous fission from Pu 244 (5). Unusual abundances of Xe 136 
have variously been attributed to extinct Pu 244 fissioning (-6) and to an 
extinct transuranium element (7). 
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In addition, some investigators detected evidence of anomalous 
alpha activity. Cherdyntsev et al. , attributed an excess of U 235 to 
either extinct Cm 147 or to an extinct transuranium element (8) , and they 
also reported an anomalous abundance of Pu 239 (9,10). They reported 
anomalous alpha activity in the 4.2-4. 6 MeV and 5. 0-5. 3 MeV regions 
for U 235 -rich samples, and in the 4. 4-4. 6 MeV region for Pu 239 -rich 
samples. 


Cherry, Richardson, and Adams (11) also reported anomalous alpha 
activity in the 4.4 MeV region, in thorite and huttonite isolated from 
Conway granite from Conway, New Hampshire. 

Gentry in his studies of pleochroic halos found "giant halos" (12), 
and an anomalous halo designated as "type Y” in a sample also containing 
Po 21 °, Po 214 , and Po 218 halos, without the halos of their usual uranium 
and thorium parents (13). 

Our investigations are being conducted by two separate techniques. 

The first technique is designed to isolate the volatile elements, because 
elements 114-118, likely candidates for stable transuranium elements, are 
analogs to the volatile elements 82-86, and would thus probably be volatile 
(3), and because the excess xenon which Anders and Heymann found was in a 
chondrite containing the volatile fraction of the elements (7). Another 
procedure’ is being used to observe non-volatile elements. 

We feel that studies of alpha activity are likely to be revealing 
as a search method. Most heavy nuclides (Z > 83) emit alphas. In particular, 
since some of our studies are on rocks containing pleochroic halos, which 
are produced by alpha particles, it is essential that we study this mode of 
emission. Alpha counting. is attractive because of its high sensitivity, 
good resolution, low background, and high counting yields for small samples. 

Our studies range from 3-15 MeV, with particular attention being 
focused on the 4.4 MeV region cited by Cherry, Richardson, and Adams, and 
on the energies corresponding to anomalous pleochroic halos, the giant halos 
representing alphas of energy 9-10 MeV (12) and the Y type of energy 6.55 MeV 
(13). In general, we first scan the entire 3-15 MeV range, and then observe 
in detail the energies less than 9 MeV, where known peaks tend to obscure 
other activity. 


II. Physical and Chemical Procedures 

A. Procedure for Non-Volatile Elements (see reference 14) 

The technique makes use of the knowledge and presumption that heavy 
elements tend to substitute themselves into the cation positions of thorite 
ETh(Si0 4 )] and zircon [Zr(Si0 4 )J . Both of these minerals are very dense, 
with a specific gravity of 6.7 for thorite and 4. 6-4. 7 for zircon (15) and 
both are non-magnetic. There is little else in igneous rock with both these 
characteristics. For example, quartz and many silicates are non-magnetic 
but have low densities. Minerals rich in iron, are heavy, but they are all 
slightly to strongly magnetic. Consequently, the technique centers around a 
heavy liquid separation followed by a magnetic separation. A flow chart is 
given below (Figure 1) . 
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The sample is first crushed in. a large manual press, such as 
the one in the Metallurgy Department of Carnegie-Mellon University, 
further crushed with a hollow tube, rod, and hammer assembly, and finally 
ground with an electrically operated mortar and pestle. Sieving is done 
at every step of the way to remove the grains which are small enough, A 
too-fine sample cannot be separated in the heavy-liquid step because of 
colloidal effects, and the grains in a too-coarse sample tend to be non- 
homogeneous. The samples used ranged from 40-200 mesh, 

A 30- g portion of a sample of the proper grain size is washed with 
distilled water, dried in an oven, and poured into a 2 50-ml’ separatory 
funnel previously filled with 70 ml of bromoform (d = 2.89 g/ cm 3 ) . The 
funnel is shaken, and the contents up to the interface with the floating 
layer are emptied into a 60-ml sintered-glass funnel over a 250-ml vacuum 
flask. The heavy and intermediate density material is dried with acetone, 
following which it undergoes an identical separation with methylene iodide 
(d = 3.325 g/ cm 3 ) . Both heavy liquids can be saved and re-used. 

The heavy fraction is further separated with a strong permanent 
hand magnet. The magnet is covered with a piece of plastic and tediously 
run over the spread-out sample as close as possible without touching. The 
attracted particles fall off the plastic onto a piece of paper when the 
magnet is removed from the plastic. 

The technique of sample deposition for alpha counting varies de- 
pending on whether or not it is desired to include volatile elements with 
the non-volatile. Both techniques use Pt discs 1 inch in diameter,, As 
much as possible of a 1-mg aliquot of the final separation fraction is dis- 
solved in 4-N HN0 3 and deposited with a dropper on the planchet. At this 
point, if the object is to count only non- volatiles , a few drops of TEG 
(tetra-ethylene glycol) are added to raise the boiling point, as the TEG 
monomer boils at 276°C and the polymers formed with heating boil even 
higher. The slow evaporation and precipitation which results gives a uni- 
form deposit (16,17). A heat lamp is placed over the planchet to dry it, 
and the planchet is then flamed to remove any remaining organic matter. 
Unfortunately, this last step, necessary with TEG, may lose the volatile 
elements. Thus when volatiles are to be included, the sample is evaporated 
from straight nitric acid. Plans are being made for setting up an electro- 
deposition- unit, so that volatile elements can be counted with better reso- 
lution than that obtained by nitric acid deposition. 

B . Procedure for Volatile Elements 


This approach is particularly useful in that uranium, thorium, and 
many of their decay products, plus most bulk elements, which tend to mask 
other alpha activity are absent from the spectrum. The sample' is ground, 
sieved, and separated with the two heavy liquids. It is then placed in 
the furnace, in a boat in a fused-quartz tube (18). The volatile material 
dissolved from the fused-quartz condensation tube with 4~N HN0 3 is ready to 
be deposited on the planchet. 
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III. Alpha Counting 

All samples were counted on one of the two ORTEC detectors available 
for this research (19). Energy calibrations were made for the narrow 
ranges with a natural uranium standard, and for the wide ranges with Sample 
D (See section V-C of this report.), which is dominated by Th 232 and its 
decay products, with lesser contributions from the U 23s series. Linearity 
within the standardized ranges was confirmed by plotting the known energies 
of the nuclides postulated to correspond to the peaks obtained, against the 
channel numbers of the peaks, and checking to see if all the points fall on 
a straight line. 

A voltage divider was added to the alpha system to extend its range. 
It was discovered that if the post amplifier bias on the alpha system is 
set at "1", the voltage divider could be used without destroying linearity. 
It was then possible to obtain energies up to 21.7 MeV. 


IV. Results 


A. Conway Granite, by Procedure for Non-Volatiles 

Samples A and B were prepared from a core sample of Conway granite, 
similar to that used by Cherry, Richardson, and Adams (11). Sample A was 
deposited with TEG and analyzed in the range 3.59-5.74 MeV. Sample B was 
deposited with plain nitric acid and analyzed in the range 3.73-5.79 MeV. 

The spectra are given below in Figures 2 and 3. Note that no activity was 
observed in the 4.4 MeV region. Sample A was run again at a wider range, 
2.00-7.56 MeV (Figure 4). 

B. Conway Granite, by Procedure for Volatiles 

The samples were obtained from a surface rock which Dr. Kohman 
collected at a quarry at Conway, New Hampshire, and which has been desig- 
nated "green" Conway granite because of its coloration. A preliminary sample 
volatilized at 800°C produced a spectrum so dominated by Fo 210 that it was 
difficult to resolve other details. It was decided to first distill out 
the Po 21 ° at a low temperature. A Sample C was prepared and heated at 600°C, 
a new quartz condensing tube was inserted, and the sample was reheated at 
700°C. The alpha spectra for Sample C at the two temperatures are shown 
in Figures 5 and 6. Po 210 was present in relatively smaller amounts in the 

second spectrum, In addition, since the overall activity of the second 
spectrum was lower than that of the first, it was suspected that most of the 
polonium volatilized at 600°C, and that the radon, radium, and bismuth are 
only partially volatilized even at 700°C. Sample C will be run at 800°C and 
900°C in the future. 

C. Chaines Anoxyennes Biotite and Monazite by Procedure for 
Non-Volatiles 


A Sample D was prepared from a biotite and monazite rock containing 
pleochroic halos, from Chaines Anoxyennes, Madagascar, obtained from Dr. Robert 
Gentry. The sample was analyzed 'at various ranges, the narrowest of which ran 
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from 0.80-8.16 MeV, and the widest from 2.1-21.7 MeV. Figure 7 shows 
Sample D analyzed' for a range of 1.085-11.505 MeV, which was narrow enough 
to preserve the details of the spectrum. 


V. Conclusions and Plans 


No anomalous alpha activity has so far been detected at this 
early stage of our studies. In the future we will look at many rock types, 
including more rocks containing pleochroic halos which Dr. Gentry sent, 
and also meteorites. 
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Figure 4. Alpha spectrum of non-volatile, non-magnetic, heavy mineral, 40-100 mesh fraction from Conway granite 
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Figure 7. Alpha spectrum of non-volatile, non-magnetic, heavy mineral, 60-140 mesh fraction from Chaines 
Anoxyennes biotite and monazite (Sample D) in range 1.085-11.505 MeV, deposited with TEG. 
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II. A. 4. 

MOSS BAUER SPECTROSCOPY OF COAL -RELATED MINERALS 
James D. Ulmer 


I. Introduction 


This project is a continuation of and expansion on earlier work 
that tried to identify iron- containing mineral constituents of coals, 
including possible "organically bound" iron compounds (1'2,3). 

The method is, to obtain Fe 57 Mossbauer spectra of various iron- 
containing compounds and minerals and to identify those in coal by 
■comparing their Mossbauer spectra. 

The parameters obtained from MUssbauer spectra that are most 
readily used for identification of chemical compounds are the isomer 
shift and the quadrupole splitting. Other parameters, line width and 
fractional absorption, for example,, can also be used to identify a 
compound. These parameters, however, depend not only on the nuclear and 
at.omic configuration of iron atoms in the compound, but also on the nature 
of the particular sample. For example, sample thickness and concentration 
of iron atoms in the sample can affect line width and fractional absorp- 
tion (4,5), and, therefore, these parameters are not so easily used for 
identification. 

In addition to the quadrupole splitting, Fe 5T displays a magnetic 
splitting in ferromagnetic materials. Observation of splitting of this 
nature, then, would help to identify such minerals or compounds. 


II. Techniques and Instrumentation 

The samples of coals, coal fractions, minerals, and compounds 
which have been utilized so far were previously prepared by Edgar Paez, 
and their preparation and chemical and thermal treatment are described 
elsewhere (2). 

The instrumentation consists of an International Chemical and 
Nuclear Corp. Model NS-1 Mbssbauer Effect Spectrometer, a Type 745 
Amperex 303 PC proportional counter, and the following Nuclear Data, Inc. 
equipment: Model ND-537 High Voltage Supply; Model ND-533 Preamplifier; 

Model ND-532 Amplifier; Model ND-522 Single-Channel Analyzer; and Model 
ND-510 Series 1100 multichannel Analyzer System used as a multiscaler. 

In addition, there are three pieces of peripheral equipment used for 
displaying and recording the Mossbauer spectra: an oscilloscope for 

immediate display of the spectrum; a Houston Instrument Model 6550 
Omnigraphic Recorder for a graphic record of the spectrum; and a 
Teletype Corp. Model 33 teletype for recording the spectrum as number of 
counts per channel on paper and paper tape. A tape reader and a card 
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punch are available in the Physics Department of Carnegie-Mellon' University 
(CMU), and they are used to convert data from paper tape to cards. 

There is also a Janis Research Company Model 6-5/8" D.T. Super 
Varitemp Dewar with Optical Tails for obtaining spectra at low tempera- 
tures, and an Eastern Scientific Company Model 4HF Laboratory Electro- 
magnet for obtaining spectra at various magnetic field strengths. 

A computer program that corrects for the parabolic shape of the 
spectrum baseline when obtained using the moving-source mode of the 
spectrometer has been written for use on the UNIVAC 1108 computer. The 
program also corrects for small, non-random differences in counting 
rates (less than 1%) for different channels and a larger difference 
(still less than 1%) in counting rates between the first 256 channels and 
the last 256 channels of the multiscaler caused by some malfunctioning of 
the multiscaler equipment. The data are corrected for these effects by 
dividing each channel of the absorption spectrum by its corresponding 
channel of a blank spectrum. Because the difference in counting rates 
between the first half and the second half of the channels of the multi- 
scaler is dependent on the counting rate, the blank spectrum and 
absorption spectrum must be obtained using the same counting rate. This 
program, then, determines the baseline and normalizes the spectrum with 
respect to this line. After the data have been corrected and normalized, 
they are printed graphically and numerically on the computer listing and 
punched on cards. 

A second program has been obtained from the Mbssbauer group in 
CMD's Physics Department that fits a Lorentzian curve to each absorption 
peak in a Mbssbauer spectrum by the method of least-squares. 


. III. Project Status and Plans 

Currently, spectra have been obtained for three samples of Jewell 
Valley coal and two samples of lignites. Spectra of these samples were 
previously obtained by Edgar Paez, and the values for the isomer shift 
and quadrupole splitting of the doublet peaks in each sample are in 
reasonably good agreement with values obtained by him. Table I lists 
the values of the isomer shift and quadrupole splitting obtained. 

These spectra have not as yet been fit by Lorentzian curves 
because the computer program obtained from the Physics Department is 
currently being revised to use our data as input. When this program is 
working, the values for the quadrupole splitting and isomer shift will be 
determined more precisely. ' 

It is believed that the peaks in the Jewell Valley coal spectra 

represent the presence of "non-pyrite iron", as designated by Lef elhocz (1). 

Further, comparison of values for the isomer shift (6) and quadrupole 

splitting (A) of the lignite samples with the values of $ and A for pyrite 
indicate that the peaks in the lignite spectra probably represent the 
presence of pyrite in the lignite samples. Values of <5 and A for the 
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coal samples investigated here can be -compared with values obtained for 
the same coal samples by Paez and with non-pyrite iron and pyrite by 
referring to Table I. 

In one of the samples of Jewell Valley coal (JVC-4 sink) 
investigated, there is a third absorption peak at approximately +0.8 mm/sec,, 
relative to the center of the sodium riitroprusside spectrum (Fig. 1). 

The fractional absorption of the third peak is approximately one-fourth 
that of the low velocity peak of the doublet. It is suspected that thi's 
third peak is:.the high-velocity peak' of a second doublet in the spectrum, 
and that the low-velocity peak of this second doublet is obscured by the 
low-velocity peak of the first doublet. The low-.velocity pealc of the 
second doublet appears to be at approximately +0.3 mm/sec., relative to 
the center of the sodium nitroprusside spectrum. In comparison, the low- 
velocity and high-velocity peaks of, pyrite are located at +0.24 mm/sec. 
and +0.84 mm/sec., respectively, relative to the center of the sodium 
nitroprusside spectrum. The second doublet observed in this sample, then, 
may or may riot be an indication of pyrite in the sample, depending on how 
well the position of the low-velocity peak was determined. Fitting 
Lorentzian curves to the spectra of this sample will determine the 
location of the low-velocity peak more precisely and it will be possible 
to decide whether the second doublet is caused by the presence of pyrite 
or not. 

Mossbauer spectra of many coal-related minerals and compounds 
have been obtained. However, none of these spectra appear to match the 
spectrum of non-pyrite iron in coal (1). This non-pyrite iron has been 
determined to be iron (II), in a high spin configuration, and, very 
probably, in octahedral symmetry (1). Mossbauer spectra will be obtained 
for coal-related minerals and compounds which may have iron with these 
characteristics. 

It has been suggested that the correlation between iron and 
silicon, aluminum, and magnesium in coal, observed by Dr. L. F. Vassamillet 
using the electron microprobe x-ray analyzer at the Mellon Institute (8) ? 
would make it profitable to obtain spectra of chlorites (2). This 
reasoning will be pursued, and various chlorite minerals will be investi- 
gated. 


The investigation of coal-related minerals and compounds at 
liquid nitrogen temperature will give additional parameters to use as 
identification because the isomer shift and the quadrupole splitting are 
temperature dependent (6,7). Furthermore, the magnitude of absorption 
can increase with decreasing temperature (5), and unobserved peaks at room 
temperature may be observed at liquid nitrogen temperature. 

It is also hoped that an external magnetic field will induce 
magnetic splitting in spectra of coals and coal fractions and thus help 
locate overlapping and unresolved absorption peaks. 
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TABLE I 


MOSSBAUER spectral parameters of coals and coal-related minerals and compounds 




This 

Work 


Other Work 

C 1 * 

Sample 

Description 

fit* 

' A™ 


6^* 


Ref. 

JVC -3 
Float. 

Jewell Valley Coal, 
physically fractionated 

1.34 + .08 

2 . 68 + 

.08 

1.41 

2.60 

(a) 

JVC -4 
Sink . 

Same, physically 
fractionated 

1 . 24 4- , 06 

2.75 + 

.06 

1.393 

2.575 

(a) 

JVC- 11 

Same, heat treatment 
200°C in. air for 3 hrs. 

1.34 + .08 

2.75 ± 

.06 ■ 

1.35 

2.575 

(a) 


Non-pyrite 
iron in coal 




1.38 

2.62 

(b) 

L-2 

Gascoyne, Grand Forks, 

N. Dak., high-iron lignite 

0.50 + .08 

0.68 + 

.05 

0.4 

0.63 

(a) 

L-3 

Beulah lignite. Hoot Lake 
Grand Forks, N, Dak. 

0.55 + .05 

0.63 + 

.04 

0.5 

0.6 

(a) 


Pyrite, Fe.S 2 




0.54 

0.60 

(b) 


# 

Designation of Paez (2). 

^"in mm/sec. 

Isomer shift ‘with respect to center of 
sodium nitroprusside. 


^Full velocity difference between paired peaks, 
(a) Paez (1967). 

^k^Lefelhocz jjt al . , (1967). 


u> 

Ln 




Relative Transmission (%) 
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4 M71- 10 SO 9 

« - * 

B. Nuclear Cosmo chemistry 


II. B. 1. 


ISOTOPIC COMPOSITION OF LEAD AND THALLIUM 
IN NATURAL MATERIALS ' 

. . James M. Huey 

I. - Introduction 

Data obtained from the isotopic analysis of lead in many natural 
materials have already found widespread interpretations in the field of 
geochronology and have been applied in the development of theories of 
nucleosynthesis. In addition, it may be possible .to extend the fine time- 
resolution of the interval following nucleosynthesis through application 
of the principles of extinct natural radioactivity to the a * 05 Pb- 205 T1 
system and isotopic analyses of thallium. Data on the isotopic composi- 
tions and concentrations of lead and thallium in stone meteorites, 
particularly the carbonaceous chondrites, should be important in the 
search for extinct natural radioactivity of 20S Pb, and lead data, when 
coupled with uranium and thorium concentration values will further the 
understanding of the history of these extra-terrestrial bodies. The 
objectives and potentialities of lead and thallium isotopic studies 
including the work of previous investigators has been described in detail 
previously (1) . 

II . Chemical Procedures 

The following scheme for isolation of lead and thallium from 
natural materials , separation-, and conversion to suitable form for mass 
spectrometry was proposed previously (1) : 

(1) Volatilization of crushed sample in gas stream to 
separate T1 and Pb from the matrix. 

(.2) Dissolution of volatiles in dilute HGl and reduction 
of Tl +3 with SO 2 . 

(3) Ion-^exchange separation of lead and thallium from 
each other and from other volatile elements . 

(4) Preparation of thallium for mass -spectrometry - as the 
nitrate. 

(5) Preparation of lead for mass-spectrometry, either as 
the nitrate or as the sulfide. 

Several additions and modifications in the chemical procedure have been 
necessary. The procedure which has been applied to the analysis of 
standard rocks,- lunar £ines, and meteorites is briefly discussed and 
outlined below: 
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S. 1 
& 


(1) Sample crushed in agate mortar and pestle to pass 
through 200-mesh nylon screen. ' 

The use of a tungsten carbide capsule and pestle 
(Wig 1 1-Bug) was abandoned when chips of the pestle 
were discovered in crushed material, thus creating 
a potential source of lead contamination. 

(2) Volatilization of sample in vacuo (~ .0.1 mm Hg) in 
fused quartz combustion tube. 

Adequate yields of lead based on analyses of _NBS 
standards AGV-1 (Andesite) and' BCR-1 (Basalt) could 
• not be obtained by heating of the crushed sample in 
a stream of oxygen. The use of a vacuum system 
produced a substantial increase in yields of both 
lead and thallium. 

(3) Dissolution of volatiles in a mixture of 5 ml 
7.9-N HC1 + 2 ml redistilled 70% HNO3 and evapora- 
tion to dryness; when concentrations are to be 
determined, spikes (92.26% Tl-203, NBS-983 lead 
std.) are added at this point. 

(4) Dissolution in T.8-M HC1 and reduction of Tl +3 

to Tl +l with S0 2 passed first through concentrated 
sulfuric acid traps for purification. 

(5) Solution placed on a 20^cm long (0.7-cm diameter) 
anion-exchange column (Dowex 1-X8, 100-200 mesh) 
and T1 eluted with 30-ml S0 2 -saturated 1.8-MHC1. 

(6) Lead eluted with 20 ml triply distilled water. . 

(7) Solutions evaporated to, dryness in quartz beaker. 

(8) Electrodeposition of lead. 

The use of a silica^el^-phosphoric-acid technique 
has been necessary for mass-analysis of small . 
quantities of lead. However, lead must be removed 
from interfering hydrocarbons by electrodeposition 
prior to mass analysis. Lead, is deposited in 
90-100% yield on a platinum anode from a 0.1-.M HNO3 
solution at 1.85 V for 24 hrs. Lead is dissolved 
from the anode in 4-N HNO3 and evaporated to dryness 
in preparation for mass analysis. 

(9) Refluxing of thallium with redistilled 70% HNO3. in ' 
quartz crucible. ' 

Hydrocarbon contamination is "reduced, and the solution 
is evaporated to dryness in preparation for mass ■ 
analysis using a triple filament-silica-gel technique. 
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. If hydrocarbon contamination in small samples, cannot 
be eliminated, refluxing with perchloric acid may -be , 
necessary.. • . 

i * 

Tests of the yields of lead and thallium obtained, incorporating 
all parts of the chemical procedure except volatilization have been made 
using small quantities of reagent solutions. The results, reported in 
Table 2, show average yields of ~ 85% for thallium and 70-75% for lead. 

It is not necessary to know yields ’’exactly since isotope dilution tech- 
niques are used, and samples are spiked prior to chemical separation, 
but high yields are desirable for adequate, recovery of the elements for 
analyses, to determine exact isotopic compositions. - ■ • <i 

These chemical procedures' have been developed in particular -for 
the recovery of the volatile elements lead and thallium from silicates,; 
metal and .troilite phases, and various separated mineral fractions of 
natural materials . At some time it may be desirable to incorporate 
these procedures into a separation scheme involving several other ele- 
ments, (e.g., Md, Sm, U, Th, Ag, Pd), which, .with the .possible exception, 
of Ag will be retained in the residue. Discussion of such a scheme is 
found elsewhere in this report (2). Some modification of the current 
procedures may he necessary for this incorporation, but basic -techniques' 
will remain the same and should be readily adaptable. 

III. Mass Spectrometry 

- Much t-ime -in this laboratory- has been devoted to experimentation 
with mass analysis of microgram quantities of lead and. nanogram quanti- 
ties of thallium. Results of these experiments are described elsewhere 
in this report (3). 

IV. Blank Analyses 

Analyses of various parts of the chemical procedure to determine 
chemical blanks for’ lead and thallium have been made at intervals 
between processing of samples. Table 1 shows the results of several 
analyses. Blanks for both' elements approach a constant value, but it 
is desirable to obtain'd, lower and more consistent lead blank. An 
analysis of various components of the blank, including reagent contamina- 
tion, quartz tubing, and atmospheric fallout contamination has been 
initiated. Blank levels for quantities of various reagents used in an 
individual separation are reported in Table I. 

V. Results 

A. MBS Standard Rocks BCR-1 and AGV-1 


T,o determine the efficiency of the entire chemical procedure 
including volatilization, several analyses of standard rocks BCR-1 and 
AGV— 1 have been made at various heating times and temperatures. In 
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Table 3, results of these analyses are reported as elemental concentrations. 
Comparison of these concentrations with reported values permits the calcu- 
lation of a chemical yield for the volatilization procedure. T1 yields are 
80% for heating times of ->_ 24 hrs . , but longer heating periods (~ 6-8 da) 
at successive temperatures are necessary to obtain comparable yields of 
lead. 

B. Lunar Fines from Apollo-11 and Apollo-12 

Preliminary results obtained in this laboratory on Apollo-11 
finds have been published recently (4). A program of differential heating 
of Apollo-11 and Apollo-12 samples' to successively higher temperatures 
was utilized to isolate various lead fractions of unique isotopic composi- 
tion. Fractions were collected at 700°, 800°, 900°, 950°, and 1100°. A 
total sample heating encompassing the entire heating range served as a 
control for comparison. Analyses of all samples have been completed. 

When detailed blank analyses have been completed, interpretations 
of the data obtained will be possible.* Lead isotopic data in conjunction 
with U and Th concentrations will aid in the determination of various 
components of the lunar fines and provide information about the thermal 
history of the moon. 

C. Meteorites 


At present, several stone meteorites are being analyzed using 
total and differential-heating techniques and the chemical procedures 
described in this report. 

The analysis of lead contents of Ladder Creek, an olivine-hyp er- 
sthene chondrite, is presented in Table 4. Mass-spectrometric data are 
corrected for a blank of 0.07 Ug. Unfortunately, a very small signal was 
obtained in the composition analysis of the high— temperature differential 
fraction (900°, 48 hr; 950°, 48 hr; 1100°, 48 hr), but a comparison 
of low-temperature (700°, 24 hr; 800°, 24 hr) and total-heating (encom- 
passing heating- range of low and high temperature fractions) samples show 
only a small compositional variation. This is apparently in contrast to 
the behavior of samples of lunar material subjected to similar treatment. 

Further examination of other meteoritic samples should prove 
invaluable in establishing- other similarities and -differences between 
meteorites and lunar samples, and indeed, those between various meteorites 
themselves. Improved lead isotopic composition data will provide infor- 
mation about the process of formation of stone meteorites and their thermal 
history. If variations in the isotopic composition of thallium can be 
observed, fine time-resolution of the interval following nucleosynthesis 
may be possible. 


Preliminary calculations indicate that variable lead isotopic composi- 
tions do exist among the various temperature fractions of both Apollo-11 
and Apollo 12. 
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TABLE I 


BLANK ANALYSES 


Procedure 

Pb Contamination 
(ug) 

Tl Contamination 
(ng) 

1 

Chemical Scheme Beginning with 

.096 

13.5 

Sample Dissolution 

i 

i 

| Chemical Scheme Beginning with 

.075 

1.48 

Sample Dissolution 

{ 

j 40 hr Heating at 950° + Chemistry 

.049 

12.02 

' Electroplating Only 

.021 

— * 

1 Heating (42 hr,' 900°; 28 hr, 975°; 

1 27 hr, 1050° C)+ Chemistry 

.084 

14.7 

! Heating (23 hr, 800°; 24 hr, 850°; 

20 hr, 900°) 4- Chemistry 

.13 

.78 

Heating (24 hr, 700°; 24 hr, 800°; 



48 hr, 900°; 48 hr, 950°; 

.17 

.64 

48 hr, 1100°) + Chemistry 



Heating (24 hr, 700°) + Chemistry 

.090 

2.58 

Heating (24 hr, 800°) + Chemistry 

.068 

.45 

Heating (24 hr, 900°) + Chemistry 

.058 

1.09 

Heating (24 hr, 950°) + Chemistry 

i 

.053 

.81 

, Heating (24 hr, 1100°) + Chemistry 
! 

.053 

.54 


Electroplating not used for Tl. 
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■ TABLE' II 

YIELDS OF COMPLETE CHEMICAL PROCEDURE FOLLOWING VOLATILIZATION 


Sample 

Electroplating Time 

Amt. Recovered (yg) 

Yield 

2.31 yg Pb ‘ 

'50 hr 

1.61 ■ 

69.7% : 

2.31 yg Pb 

50 hr 

1.67 

72.3% 

2. .31 Ug Pb 

24 hr 

1..68 

72.8% 

2.31 Ug Pb 

24 hr 

2.17 

93.9% 

6.04 Ug Tl 

None 

5.23. 

86.6% 

6.04 yg Tl 

None 

5.37 

89.0% 





TABLE III 


•P- 

-P- 


MBS STANDARD ROCK ANALYSES 


Sample 

Amt. (g) 


Treatment of Sample to Remove 
Volatiles Prior to Chemistry 

Amt . Recovered 

(hg) 

Pb Tl' 

Derived Mass 
Fraction 
(ppm) 

Pb Tl 

Yield 
Pb ■ 

<XS 

Tl 

BCR-1* 

2.77 

Heat 

(1200°, 3 hr in 0 2 ) 

2.02 

.230 

.73 

.083 

5.4 

25.1 

BCR-1 

2.65 

Heat 

(1100°, 4-1/2 hr in Oj 

- 2.05 

.347 

.77 

.131 

5.7 

39.8 

BCR-1 

1.39 

Heat 

(1000° for 3 hr; in 0^ used cold trap) 

4.02 

.503 

2.90 

.364 

21.0 

110.0 

BCR-1 

2.78 

Heat 

(1000° for 3 hr; in 0 used cold trap) 

3.44 

.344 

1.24 

.124 

9.2 

37.6 

BCR-1 

2.47 

Heat 

(1000° for 3 hr; in 0^ used cold trap) 

.66 

.329 

.27 

.134 

2.0 

40.5 

BCR-1 

3.023 

Heat 

(1000° for 3-1/2 hr; vacuum) 

14.7 

.640 

4.86 

.212 

35.8 

64.4 

BCR-1 

1.995 ' 

Heat 

(1000° for 4 hr; vacuum) 

21.2 

.422 

10.6 

.212 

78.3 

64.2 

BCR-1 

1.082 

Heat 

(900° , 10 hr) 

11.7 

.282 

10.8 

.260 

79.6 

79.0 

BCR-1 

3.61 g 

Heat 

(900°, 30 hr) 

42.3 

' .754 

11.7 

.209 

86.5 

63.4 

AGV-1^ 

1.817 

Heat 

(1000°, 24 hr) 

18.8 

.534 

10.4 

..294 

28.4 

86.1^ 

AGV-1 

2.000 

Heat 

(900°, 72 hr) 

43.2 

.622 

21.6 

.311 

59.2 

94.5 

AGV-1 

2.502 

Heat 

(800? 20 hr; 900? 20 hr; 1050? 20 hr) 

43.1 

.750 

17.2 

.299 

47.1 

87.5 

AGV-1 

. 2.630 

Heat 

(900°, 24 hr) 

23.7 

.750 

9.0 

.285 

24.7 

83.5 

AGV-1 

2.009 

Heat 

(800? 20 hr; 900? 20 hr; 1050? 20 hr) 

65.8 

,645 

32.8 

.321 

89.8 

94.0 

AGV-1 

3.227 

Heat 

(800? 20 hr; 900? 20 hr;, 1050? 20 hr) 

31.8 

.905 

9.9 

.280 

27.0 

81.9 


Yields of BCR-1 based on 13.56 ppm Pb(6) and .330 ppm Tl(7). 


j. 

r Yields of AGV-1 based on 36.53 ppm Pb(6) and .342 ppm Tl(7). 



TABLE IV 

LEAD ISOTOPIC DATA OBTAINED FROM DIFFERENTIAL 
HEATING OF LADDER CREEK METEORITE 


Sample 

' Mass (g) 

206/204 

207/204 

208/204 

207/206 

208/206 

Pb Mass Fraction 
(ppm) 

Total Heating 

(Sum of Low + High Temp.) 

3.994 

19.63 

16.30 

38.74 

.8304 

1.9733 

.518 

Low Temp. Fraction 
(700°, 24 hr; 800°, 24 hr) 

5.300 

19.08 

15.96 

38.47 

.8363 

2.0158 

.123 

High Temp. Fraction 
(900°, 48 hr; 950°, 48 hr; 
1100°, 48 hr) 

5.300 






.226 

$ r- 
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II. B. 2. 


CHRONOLOGY OF GALACTIC HEAVY-ELEMENT NUCLEOSYNTHESIS 


James M. Huey and Truman P. Kohman 


I. Introduction 


Three quite independent methods for obtaining quantitative esti- 
mates or lower limits of the n age of the universe" have been discussed 
by Fowler (1) and Sandage (2) , with the conclusion that each of these 
methods yields results consistent with an age between about 8 and 25 Gy: 

(a) calculation of model ages of the oldest stars from the main-sequence 
cut-off points in color-luminosity diagrams for star clusters; (b) deter- 
mination of model parameters for a monotonically expanding universe from 
galactic red-shifts; and (c) evaluation of parameters in models of galactic 
nucleosynthesis from relative abundances of radioactive and stable products. 

We have undertaken a refinement of the third method, using chiefly 
a model for galactic nucleosynthesis proposed earlier by one of us (3,4) 
with recent values of cosmochemical and nuclear -physical data. 


II. Multi-parameter Models of Galactic Nucleosynthesis 

In models of the time-dependence of nucleosynthesis, the number 
of astronomical parameters which can be evaluated is equal to the number 
of critical observable quantities which can be related' to them through 
certain nuclear parameters. A number of models having different' numbers 
of astronomical parameters have been proposed, and these are summarized 
in Table I. The numerical part of each Model Number is the number of 
astronomical parameters. 

The symbolism is an extension of that used earlier (3) : 


Capital Greek letters : 
Small Latin letters: 
Small Greek letters: 
Capital Latin letters; 


Time intervals . 

Specific times . 

Model and nuclear parameters . 

Extensive quantities related to matter. 


The symbols in Table I and the ensuing discussion have the following 
significance: 

t = time 

o - present time (t = 0) 

m = time of meteorite formation (as cool solids capable of 
retaining rare gases) 
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s = time of beginning of solar-system formation (defined as end 
of galactic nucleosynthesis of any solar-system 
matter-, presumably the time of eruption •of the 
last contributing r-process. Type I> supernova) 

f = time of beginning' of. formation of r-process nuclides 

g = origin of our galaxy (and of the universe) 

M = o - m = age of meteorites (as cool solids) 

S = m - s = decay interval for extinct natural radioactivity 

(especially with radiogenic Xe) 

E = o - s = age of solar system = E + M 

.A = s - f = duration of galactic nucleosynthesis of heavy 

elements of the solar system 

'E - *f g = average evolutionary time for Type I supernovae 

r = o — g = age of galaxy = E + A + d+ M 

A = mass number of nuclide 

N. = number of atoms of nuclide of mass number A 
A . 

. =- cumulative r— process production rate of nuclide of mass 
■ number A 

A^ = decay constant for radionuclide of mass number A 

A = decay constant for general galactic stellar activity 

:s (assumed proportional to r-process nucleo synthetic 

activity) 

cp = fraction of total' r-process elements in terminal spike 
(from last Type I supernova contributing 
to solar system) 

a = fraction of total r-process elements in -initial spike 

(from massive" first-generation .stars during collapse 
of galaxy to disc configuration) 

Xjj = (N2 3 5/N2 3 8') 0 /(P2 3 5/P2 3£t) 

X T Jj ~ (N232/N238) 0 /(p232/P238) 

Xj = (Ni29/Nl27) m /(P129/Pl27) 

Xp u = (N24if) m /(N23 8) 0 /(p2«tfr/p2'3 8) 
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As explained earlier (3), the abundance ratio of two isotopes or otherwise 
related nuclides always enters into the mathematical equations as the ratio 
to the corresponding production-rate ratio, and we have adopted Clayton's 
(5,6) designation of such ratio-ratios by the symbol x with a subscript 
for the element in the numerator of each ratio. 


III. Evaluation of Parameters 


The general equation for each pair of related nuclides is 


'■El 


( W<V„ 


Here we have adopted the philosophy of relating the astronomical parameters 
of a nucleosynthesis model as directly as possible to observed quantities . 
Thus, the abundances of stable and extant radioactive nuclides are taken 
at the present time; for extinct natural radionuclides the time x in the 
above equation is taken as m, which is when the observable abundances of 
the radiogenic products (Xe isotopes in all cases so far) were established. 
X - o - x is the corresponding interval, 0 or M. 


In addition to the four critical observable isotope or nuclide 
ratios already mentioned, there is a fifth, which may formally be taken 
as the ratio of radiogenic Pb 207 to radiogenic Pb 206 produced from 
uranium in meteorites during the interval M. Actually, there is a large 
body of evidence based on five radioactive systems (U 238 -Pb 206 , U 235 -Pb 207 , 
Th 232 -Pb 20 8 , Pb 87 -Sr 87 , K 40 -Ar 40 ) which give more-or-less concordant 
estimates of M, the age of the meteorites (the most accurate of which 
actually does come from radiogenic Pb 207 and Pb 206 via the slope of the 
meteoritic Pb 207 /Pb 204 versus Pb 206 /Pb 204 isochron). 


One could set up a fifth equation for the radiogenic Pb Z07 /Pb 206 
ratio and solve the five simultaneous equations for the five parameters 
of Model 5a or 5b. However, the solution of the equation for M is com- 
pletely independent of the other equations and input quantities, so we 
choose instead to accept the previous estimates of M by other workers, 
and consider M instead as an input quantity. This leaves four simul- 
taneous equations (for Xrj> X Th > Xt> X p ) involving four additional 
parameters to be solved for (see x Table I) . 


As the "nominal" value of M, we have taken 4.56 Gy, which is the 
best value of the age of chondrite meteorites from Pb-Pb isochrons (7). 
The nominal values of the other input quantities , with the abundance 
ratios and production-rate ratios from which they were derived, are 
listed in Table II. According to recent observations of Podosek (8), 
there are definite differences in the ratio of radiogenic Xe 123 to stable 
I 12 in different chondrite and achondrite meteorites, indicating dif- 
ferences in formation (cooling) times amounting to as much as 50 My. 
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Therefore, for consistency, we have used values of ( 1 ^ 29 /^ 127 ) and 
(N 244 )^/ (N 2 3 8 ) ^ derived from the same meteorite, St. Severin (whole- 
rock) in the same laboratory (8,9). 

The four simultaneous equations have been solved using the CMU 
UNIVAC 1108 computer and a FORTRAN V program designated CHRONS. The 
program provides solutions for either Model 5a or 5b. It is based on 
the Newton-Raphson method of successive approximations starting with 
estimates of the values of the parameters. Each input quantity, in- 
cluding M, is systematically varied through the physically interesting 
or possible range on both sides of its "nominal" value, and separate 
solutions are obtained. Variations of a X correspond to variations in 
either the abundance ratio or the production rate or both. 


IV . Results and Discussion 


The "nominal" values of the parameters for the two models are 
given in Table III. Figures 1 and 2 show how the output quantities vary 
with the input quantities for Models 5a and 5b, respectively. In each 
plot, each of the curves shows how two of the output quantities vary with 
one of the input quantities. The numbers on the X curves represent the 
ratio of the input value to the nominal value. The numbers on the M curve 
represent the absolute values used for the age of the meteorites expressed 
in Gy. 


We feel that Model 5a is the more realistic of the two 5-para- 
meter models which have been proposed so far. Contrary to some opinions, 
we feel that this model, with parameters in the ranges of derived values, 
is in satisfactory concordance with astrophysical data on the build-up 
of heavy elements in the galaxy, with the ages of galactic halo and disc 
star clusters, and estimates of the age of the universe based on galactic 
red-shifts. Unfortunately, the derived value of A, the duration of 
galactic nucleosynthesis, is very sensitive to and hence to both 

the Th/U ratio and P 232 /P 238 * Therefore, precise statements about the 
age of the universe (F) by this method cannot be made. Similarly, the 
value of A., or a depends strongly on both Xj^ and X-y* Thus, to obtain 
information about the early phases of galactic nucleosynthesis will require 
very good observational and theoretical information about the very long- 
lived radionuclides, U 235 , U 238 , and Th 232 , whose present abundances re- 
flect in part the early time-dependence of r-process nucleosynthesis in 
the galaxy. 

After completing these calculations, we learned that Fowler (10) has 
made an almost identical set of calculations for what we have designated 
as Model 5a. His "standard" input values differ somewhat from our 
"nominal" values, so that his best output values differ somewhat from 
ours. His ranges of values for the parameters overlap ours, and confirm 
that A and A }{ are very uncertainly determined by this method. His 
"standard" results cp = 8.8% and E = 175 My are both considerably higher 

than our "nominal" values, and we are- investigating the reason for this 
apparent disagreement. 
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We feel that neither Fowler's "standard" input quantities nor 
our "nominal" values .are the best that can be used at the present time, 
and we are continuing with an evaluation of all available abundance data 
and production-rate estimates as a preliminary to repeating the calcula- 
tions for both Models 5a and 5b. 
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Table X 


MODELS FOR THE TIME -DEPENDENCE OF GALACTIC NUCLEOSYNTHESIS OF SOLAR-SYSTEM HEAVY ELEMENTS 

Ui 

NO 


Model 

No. 

Description of Model 

Reference to Model 

Parameters 

Input Quantities 

1 

Single -event 

Houtermans 1947 ( 11 ) 

2 



*U 



2a 

Continuous uniform 

Burbidge, et_ al . 1957 (12) 

2 

A 





2b 

Single -event 

Reynolds i960 (13) 

M 

M 
— i 


[M] Xj 



3 

Continuous uniform 

Wasserburg, et_ al . i960 (l 4 ) 

M 


A 

EM) Xj 


(* * * ) 

4 a 

Continuous decaying 

Fowler and Hoyle i960 ( 15 ) 

M 

W 

D> 

P' 

(M] Xj X„ X Ih 

[MHXj] X u 0 s l87 /Re 187 

(t) 

4 b 

Continuous uniform with 
initial spike 

Dicke 1962 ( 16 ) 

M 

H 

M 

A a 

[M] [Xj] 


(ft) 

4 c 

Continuous uniform with 
terminal spike 

Wasserburg, et_ al_. I969 (17) 

M 

H 

A 9 

[H] Xj. X pu X„ 

X Th 

(§) 

4 d 

Initial and terminal spikes 
with no continuous 

Wasserburg, et_ al. 1 969 (17) 

M 

77 

A 9 

'“1 X I X Pu X U 

X Ih 

( § ) 

_ 5 a 

Continuous decaying with 
terminal spike 

Kohman I96I ( 3 ; 4 ) 

M 

H 

A X.. 9 

[“) X ! X P U *0 

X Th 


5 b 

Continuous uniform with 

initial and terminal spike 

Hohenberg 1969 (l 8 ) 

M 

M 

A a 9 

X I X P U \ 

X Th 

(§) 


* Model not completely evaluated. 

f- Clayton 1964 ( 5 > 6 )j system over -determined and consistent solution not obtained, 

ft Model not independently evaluated. 

§ Unique solution not obtained. 




Table IX 

NOMINAL INPUT QUANTITIES FOR SIMULTANEOl 


Abundance Ratio 

^235/^238^° - 

.00726 

11 

0 

^~CO 

OJ 

23 ^ 

3-66 

^129^127^ 

. 0000775 

( N 244 V^ N 238 ^o 

= .0256 

rad 


( N 2oy/ N 2 06 ) m _ 0 

= . 3946 


Production-Rate Ratio 


P235/P238 = 1 *^ 
^252^238 = 1 * 7l| ‘ 
Pl29 /l ^127 = 1 
p 244 //p 238 = 


Ref . 

Resultant . 

( 19 ) 

x„ ■ = .00^97. 

. ( 19 ) 

*!h = 2 ' 10 


X = .0000773 

( 19 ) 

x Pu ' -°5 8 . 


M =' 4.56 Gy 
















Table III 

NOMINAL VALUES FOR ASTRONOMICAL PARAMETERS 


Quantity 

Model 5a 
(Kohman I 96 I) 

Model 5b 
(Hohenberg I 969 ) 

[M] 

4 . 56 Gy 

4,56 Gy 


£ 

-=f 
« — 1 

142 My 

<P 

3*5$ 

2.c$ 

a 

— 

65$ 


0-6l Gy 1 

— 

A 

3.92 Gy 

3-17 Gy 
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Figure 1. Astronpmical parameters and their dependence on input 
quantities for galactic nucleosynthesis Model Ho. 5a 
(continuous decaying nucleosynthesis with terminal 
spike) . 
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Figure 2. Astronomical parameters and their dependence on input 
quantities -for galactic nucleosynthesis Model No . 5b 
Ccontinuous uniform nucleosynthesis with initial and 
terminal s-niV^cS . 
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C. Nuclear Chemistry 


II. C. 1. 

search for natural radioactivity OF CALCIUM- 4 : 8 

Mark J. Yeager 


I. Introduction 


In searching for single beta decay of Ca 48 , we are employing the 
active daughter extraction technique proposed by Jones and Kohman (1,21: 
Sc 48 activity is chemically isolated from large amounts of calcium. Thus 
far the work has been concerned with - the development of a procedure which . 
removes any radioactive impurities and which can be conducted in a time 
which allows detection of the 44-hr half-life of Sc 48 . 

Rajagopalan (3) has recently placed a lower limit of 0.9 X 10 2 ° yr 
for single beta decay of Ca 48 . This value is an increase by a factor of 
36 over the latest limit determined in this laboratory by Van Itallie 
et al. (4), 2.5 x 10 ls yr. In his investigation, Rajagopalan employed a 
g-y coincidence counting system, and the Sc 48 activity was coprecipitated 
with Fe(OH) 3 from ~ 240 kg CaCl 2 (78 kg natural Ga) . 

In order to improve Rajagopalan’ s estimate or detect Sc 48 activity, 
counting sensitivity and the amount of calcium' used are critical.. The 
large-scale experiments will be conducted using ~ 100 kg of natural calcium 
(184 g Ca 48 ) . A Sharp Laboratories Model LB-100 LowBeta proportional 
counter with a background of .8 + .2 cpm has been utilized for all sample 
counting. The Sharp system contains two side-window Geiger .counters, 
each mounted directly below a guard Geiger counter. The new detectors 
installed in 1965 have a window diameter of 2-1/4 inches (5) . Samples 
are counted on aluminum planchets with a 5-cm diameter. 

Five 1-Curie plutonium-beryllium neutron sources have been used . 
to promote neutron-capture reactions of Ga 46 and Ca 4s . The Sc 47 activity _ 
(3.4-day half-life) and the Sc 49 activity (1-hr half-life) serve as 
radioactive tracers for testing radiochemical procedures (6) . The 
reactions involved are (7): 


Ca 46 (n,y) Ca 47 
Ca 48 (it, y) Ca 49 


e > 

t-L /2 = 4.53 day 

e: > 

1 1 / 2 = 8.70 min 


Sc 47 


Sc 49 


E > 

t l/2 = 3.35 day 

fil > 

t i / 2 = 37.5 min 


Ti 47 


Ti 49 


The procedure for homogeneous irradiation has been discussed by Fang (6) . 



II. Dependence of Scavenging Efficiency of Sc Activity on pH 

Several initial experiments have been conducted using the neutron 
sources to determine the optimum pH for Fe(0H) 3 scavenging of scandium 
activity. 4 £ of a 5-M Ca(N0 3 ) 2 solution were irradiated for one hour, 
after which .125 g Fe' 1 "*' 1 ' was added. The Fe(OH) 3 precipitate was cen- 
trifuged, and a second precipitation from the Ca(N0 3 ) 2 solution was made 
using '.125 g Fe" H_f \ The precipitates were dried on separate planchets 
and counted. The time for extraction of both precipitates was ~ 4 hours. 
A graph of initial activity of Sc 49 versus pH indicates that the optimum 
pH for scavenging is 7.5 + .1 (Fig. 1). At pH < 7.5 the Sc(OH) 3 has 
probably not completely precipitated; it has been reported that Sc(OH) 3 
does not begin to precipitate until pH 7 (8). And at pH > 7.5 Sc ++ "*‘ has 
probably complexed with ammonia,/ which was used to adjust the pH of the 
Ca(N0 3 ) 2 solution; the existence of a hexamino scandium cation in the 
presence of excess ammonium ion has been established (8). Because of .the 
possibility of such complexing, Ca(OH) 2 has been used to raise the pH in 
all succeeding experiments. 


III. Scavenging Efficiency of Fe(OH) 3 

To determine the scavenging efficiency of Fe(OH) 3 , three successive 
precipitations were made using a 4 i aliquot of 5-M Ca(N0 3 ) 2 solution (pH 
7.5) and .125 g Fe+++. The initial Sc 49 activities for the three samples 
were respectively 639 + 10 .cpra, 70+7 cpm, and 12 + 6 cpm. Assuming that 
the sum of these three activities represents virtually 100 percent of all 
Sc 49, activity, the recovery (from the first precipitation is 88.8 + 2.2%. 


IV. Radiochemical Purification Procedure 


The main structure of the procedure was provided by Van Itallie 
e_t al. (4,9) and Fang (6), and their procedures originated from The 
Radiochemistry of the Rare Earths, Scandium, Yttrium, and Actinium (8) . 
The procedure was perfected using 4-i aliquots of 5-M Ca(N0 3 ) 2 solution 
(pH 7.5) irradiated to saturation activity with the neutron sources. 

.125 g Fe" ,-H " was used to precipitate Fe(0H) 3 , and the solid was collected 
by centrifuging and washed with ~ 10 ml H 2 0. After dissolving the iron 
in 9-M HC1, 10 mg scandium carrier was added, and the Fe' was removed 
as FeCl 3 using a 15-cm anion-exchange column with a diameter of 2.2 cm 
(Dowex 1-XtO ; 200-400 mesh) . The effluent HCl was evaporated to ~ 10 ml 
and diluted with H 2 0 to ~ 50 ml. Drops of 6-M NaOH were then added to 
precipitate Sc(0H) 3 . Because the sodium salt of ScF 6 is insoluble in 
aqueous solution, the solid was washed twice with ~ 5 ml to remove 

'any traces of Na + . After dissolving the Sc'(0H) 3 in drops of 6-M HN0 3 , 
the volume of the solution was doubled with H 2 0, and an equal volume of 
6-M NH 4 F was added. By adding 10 mg lanthanum carrier and heating the 
solution, LaF 3 and any rare earth contaminants such as YF S precipitated, 
and the scandium remained in solution as ScF 6 . The volume of the 
solution was then doubled with concentrated HC10 4 and heated to keep the 
majority of the insoluble ammonium perchlorate dissolved. ScF 3 , the 
resulting precipitate, was isolated by centrifuging and washed with 
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~ 20 ml H 2 0 to remove any traces of ammonium perchlorate. The ScF 3 was 
dissolved using ~ 5 ml saturated boric acid and ~ 2 ml concentrated HN0 3 . 
Drops of 6rM NaOH were added to precipitate Sc(0H) s , and the solid was 
washed twice with ~ 5 ml H^O and once with the same amount. of acetone. 

The precipitate was transfered to a counting planchet using acetone to 
facilitate drying under an infrared lamp; the sample was then counted. 

The average time for purification was ~ 5 hours. 


V. Results and Conclusions 

To indicate the success of the procedure, unirradiated experiments 
were conducted, and the blank samples were counted. Some short-lived 
activity and a long-lived component with an activity of .2 cpm are 
indicated (Fig. 2). Possibly these impurities will be removed after 
several successive Fe(0H) 3 precipitations. In order to improve the 
procedure, a TTA extraction is presently being added as an additional 
purification step, and triply-distilled B 2 0 will be used in future 
experiments. 

Both 57.5-min Sc 49 activity and 3.4-day Sc 4 " 7 activity were 
detected in the irradiated samples for testing the chemical procedure • 
(Fig. 3). Since both of these scandium activities are generated by 
beta_ decay of calcium, they would be in the same chemical state as Sc 48 . 
•Thus, it is reasonable to consider the chemical procedure applicable for 
Sc 48 extraction and purification. By comparing the initial activities 
for Sc 4 ^ between a purification experiment and for the experiment in 
which three successive Fe(0H) 3 precipitations were made, a value of 
~ 16% is obtained for the overall chemical .yield. 


VI. Future Large-Scale Experiments 

The Ca(N0 3 ) 2 solution used by Van Itallie at ^l._(4), has been 
homogenized and purified by Fe(0H) 3 precipitation. Because NH 4 0H was 
used to adjust the pH, it is possible that a significant portion of any 
scandium present would be complexhd. For this reason a fresh solution of 
Ca(N0 3 ) 2 is being prepared. 

Because of the. high Ca(N0 3 ) 2 concentration in Van Itallie' s 
solution (5-M), Fe(0H) 3 settling is severely inhibited: instead of 

taking a few hours for the precipitate to settle in the 50-gal tanks, it 
takes ~ 5 days. A Helixtractor manufactured by International Equipment 
Company was fitted to the centrifuge, but the device failed to collect 
the precipitate. An alternative procedure is presently being investi- 
gated: larger amounts of Fe’hH' will be used to scavenge the scandium 

from a more dilute solution of Ca(N0 3 ) 2 ; and after collecting the 
precipitate, the Fe'~^ will be removed by anion exchange. If anion 
exchange fails to remove all the iron, the majority of the Fe+++ will be 
separated by solvent extraction with di -isopropyl ether. The remainder 
of the iron will then be removed by anion exchange. In recent experiments 
the precipitation time has been reduced to 20 hours by using 1.0. g F6+++ 
per 50 gal 4-M Ca(N0 3 ) 2 solution. To coagulate the precipitate quicker, 
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the solution containing the Fe 1 1 •' [Fe(N0 3 ) 2 • 4H 2 0 in 1-M HN0 3 ] has been 
heated to 70°C before addition to the Ca(NO s ) 2 solution. Immersion 
heaters will eventually be tested as a means of reducing the settling 
time by inducing coagulation of the Fe(0H) 3 . 

i 

A single extraction experiment using Van Itallie's solution has 
been conducted, but no 44-hr Sc 48 activity was detected (Fig. 4) . After 
determining the proper conditions- for Fe(0H) 3 settling, large-scale 
experiments will be carried out in search of single beta decay of Ca 48 . 
More elaborate counting conditions will eventually be employed to improve 
sensitivity. 
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Figure 1. Dependence of scavenging efficiency of scandium 
activity on pH of 5-M Ca(H03 ) 2 solution. 
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Figure 3* Decay curve of purified Sc( 0H)3 extracted from k liters of 
5-M Ca(NCb) 2 solution irradiated to saturation activity 
with 5-Curie Pu + Be neutron source* 
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Time after Isolation from 5-M Ca(N03 ) 2 Solution (days) 

Figure h. Decay curve of purified Se(QH) 3 extracted from 
~ 80',kg calcium. 
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INSTRUMENTATION 
Mark Haramie 
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A brief description of all equipment and plans for modifications and 
acquisitions will be made. Only equipment not mentioned in previous reports 
and new modifications will be described in detail. For further information 
concerning equipment as of the beginning of the year, refer to Toia (l). 


A. Mass Spectrometry 

A Nuclide Corporation mass spectrometer Model SU 2.4 (2,3) is available to 
measure abundances of isotopes and their ratios. The instrument has a triple - 
filament solid -source, with an analyzer tube of a 12-inch radius and 90° deflection. 
A resolution of less than 1 mass unit at mass 300 can be obtained. Sensitivity 
can be increased with a 20-stage "ion-electron multiplier" collector. 


1. AC Power; Upon original installation, regulated power was supplied to 
the magnet -accelerator assembly, vibrating -reed electrometer, and ion -electron 
multiplier to reduce short-term drift. The green power strips located inside the 
back of the instrument are the regulated power lines. The gray are the unregulated 
power. 


2. Magnetic field : The instrument was supplied with a magnetic field that 

can be changed with either a sweep control or stepping toggle switch. The sweep 
control changes the magnetic current at a preselected constant rate; or manual 
operation of magnet current can be used. 

The stepping operation was a modification on the original system. Previously 
the magnet current was controlled by the difference in voltage between the center 
taps on two potentiometers in series of 20-k £2 and 200-2 which are also in series 
with a 1300-2 resistor. To this was added a l-k2 potentiometer in series with the 
voltage -divider network shown in modified electronic sweep control for magnet 
current diagram. Figure 1. The additional resistance, determined by the setting of 
the potentiometer, will change the magnet current ten times- as much as the 200-2 
potentiometer. There is a switch which will bypass the new resistance or leave it 
in the circuit, enabling a change in magnet current instantly. 

A limiting factor of slow response in changing magnetic field arises because 
of hysteresis losses. When a magnetizing force is applied to a ferromagnetic 
material there Is always a lag in magnetic intensity. The time to reach equilibrium 
depends on the amount of change in the magnetizing force, in this case the current. 

A waiting period of four to forty seconds is needed to change from one mass to another. 

Solution to the problem of slow response times for magnet equilibrium is to 
provide a detector to sense the magnetic field. The detector would then send a 
signal back to a unit that would compare a reference with the actual signal and then 
make a correction in magnet current. A probe field detector with feedback is the 
type of device which will do the operation. Plans are being considered to 
incorporate a detector of this type in our present system. 
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3. Digital Readout : During the past year a digital -output system from 

the vibrating-reed electrometer was acquired and installed. It consists of a 
Hewlett Packard 2212A voltag| -to-frequency converter with a zero stability of 

maximum nuiiiber^of -150,000 pulses full-scale is achieved with a linearity 
Coirbpui pluses 3 ft dn the converter are fed to a Monsanto 1510 A counter -timer. 
Dedd-t'ime of the counter is not important because the pulses from the converter are 
evenly spaced. The counter is equipped with a ti m er that will open the counter gate 
and allow pulses to be counted for periods of .1, 1, or 10 seconds . Display of the 
digital format is either by lit numbers on the face of the counter timer or by a 
Hewlett Packard 5050 B digital recorder. Numbers from the recorder are printed in 
sequence on a continuous paper strip. 

A paper -tape punching unit or a computer -card punch is being considered. This 
will allow direct submittal of data to a computer. Already a computer program is 
being used to calculate isotopic ratios, so the addition of such a unit will greatly 
increase the efficiency of the projects involved with the use of the mass spectrometer. 

B. Radioactivity^ easurement Instruments 

1. Ge(Li) Gamma -Spectrometry System : The Germanium -Lithium gamma -ray 

spectrometer system (1) has been rimningTthe past year with only a repair on the 
512 -channel analyzer. During this time a l/8-inch copper shield was added to the 
inner walls of the lead shield to reduce the effect of 8l-keV X-rays coming from 

the shield. The X-rays originate from the effect of lead fluorescence resulting from 
cosmic radiation, natural radioactivity, and sample radiations. The addition of the 
copper reduced the activity low enough to not he detected in an 80 -minute count. 

Also, a plexiglass frame was constructed to support samples at different distances 
above the detector. Sample placement at 6 inch intervals with five different 
locations enabling a maximum distance of 3 feet is obtained. 

2. Alpha- Spectrometer Systems : Alpha counting to the present has been done 

on two identical systems (1) using 200 -mm 2 0RTEC detectors. Recently two ^50-ram 2 

silicon surface barrier detectors have been purchased and installed. The present 
system resolution of system number 2 is. 32 keV FWHM. System number 1 seems to have 
a problem located in the preamplification, causing resolution to fall off to 6o keV 
FWHM. Maintenance has been started to correct this. 

Amplifiers provided with the system enable a maximum width of energy of only 
5.5 MeV to be seen at one time using our 256 -channel analyzer. By the use of a 10-k£2 
potentiometer in parallel with the signal from the amplifier, a smaller signal was 
taken from the amplifier to increase the energy range that could he seen at one 
time. For a discussion of the energy ranges obtained refer to the section by Landes 
in this report (4). 

C. Mossbauer Spectrometry 

In last year's report (l) it was stated that an International Chemical 
and Nuclear Corporation Mossbauer spectrometer Model MS-1 was purchased by the 
Chemistry Department, hut there seemed to he a problem with unequal dwell times 
in the different halves of the memory associated with the Nuclear Data Model 1100 
Multichannel Analyzer. Random or pulse -generator pulses fed directly into the 
multichannel analyzer or fed into the beginning of the electronic Mossbauer network 
produced a non -random baseline that has repetition with periodicities. Also the 
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average number of counts in the second half' of the memory has a greater amount 
then the first half producing a step at the 256-257 channel. Originally the step 
was greater than 1$ but has been reduced to ~ 0.2 $ by Nuclear Data. This is still 
not satisfactory and arrangement s are now being made to have the unit or complete 
system replaced in the near future. 


In the meantime, a computer program has been written to divide each spectrum 
with a separately determined baseline. Data collected and determined by this method 
seems to be satisfactory for spectra with relatively few and strong absorption peaks(5)« 


Unmodified Systems and Equipment 

The following equipment and systems are available to project personnel. 
For further description of equipment and systems refer to reference (l). 

Sharp Low -Beta Counter Model LB 101 

Gas -flow proportional counter 

Harshaw Model 858 "2 x 2" Hal(Tl) detector 

4 -pi well -type "k x 4” Nal(Tl) detector 

Coincidence system employing two Harshaw detectors Model 12SC12-H "3 X 3" Nal(Tl) 
detectors . 

Coincidence system employing two Harshaw 8s8 u 2 X 2” Nal(Tl) Integral-line 
scintillator detectors . 
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GENERAL SCHEME OF CHEMICAL 'ANALYSIS OF LUNAR AND METEORITIC SAMPLES 


H. Ihochi 


I. Introduction 


Basing on the development of space research, especially with 
respect to theories about the origin of the solar system and cosmic 
abundances of the elements,, research on extraterrestrial materials has 
gained considerable importance. Because of limitations in available 
quantities of such samples, however, a development of a systematic analyt- 
ical scheme for the determination of the elemental abundances and isotopic 
compositions of a considerable number of elements from a single sample .of 
a complex matrix is required. . As a part of the continuing investigation 
at Carnegie-Mellon University of the possible existence of additional ex- 
tinct natural radioactivities, such a general scheme has previously " been 
proposed by J. .Tanner (1) for the isolation and determination of the ele- 
ments, Pd, Ag, Pb, Tl, Sm, Nd, U and Th in lunar materials. Possible 
extension of this to include W, Hf, Ca, K, Mg, and A1 has been discussed 
by Kohman (2) . 

I-n the following section, recent modifications and application 
of this scheme are discussed. 


II. Experimental and Results 

A. Separation of Lead and Thallium by Volatilization 

It has been demonstrated that the previously proposed isolation 
method for separation of Pb and Tl by volatilization in an oxygen flow 
system as proposed by Huey (3) yielded rather low chemical yields for both 
lead and thallium. By means of a vacuum volatilization, considerable in- 
crease- in chemical yields were obtained, and oxidation of sample was avoided. 
The diagram of the vacuum volatilization system is shown in Figure 1. The 
results obtained. in the above experiments are discussed in detail elsewhere 
in this report (4) . 

.The behaviors .of silver and rare earth elements in this volatili- 
zation step were tested by radioactive tracers Ag 110 and Sm 145 . The results 
indicated that silver was partially volatilized and that rare-earth elements 
were quantitatively retained in the non-volatile residue. The fraction of 
silver activity in the volatile condensate was found to be very dependent 
on the condition of sample. When only the silver tracer was used, approxi- 
mately 100% of the silver activity was found in the volatile fraction. When 
the same tracer was mixed with basaltic rock, the fraction, in the volatile 
deposit varied from 30% to 70% depending on the operating temperatures of 
800°C to 1000 °C , respectively. 



* Decomposition and Dissolution of Non-Volatile Residue 


Complete decomposition of the non-volatile residue was attained 
when the residue was repeatedly treated with a mixture of concentrated 
HF, concentrated HCl, and concentrated HN0 3 of approximately equal volumes. 
This decomposition procedure was carried out in a teflon beaker covered 
with a teflon plate at temperature slightly below the boiling point of the 
mixed acid. Extreme care must be taken in order to avoid complete dryness 
of the sample during this treatment. Complete dryness of the sample tended 
to produce insoluble residues. 

Treatment of the non-volatile residue with HF- HN0 3 mixture or 
HF-HCIO4 mixture resulted in formation of insoluble oxides, most probably 
aluminum oxide. Insoluble sulfate was found when the decomposition was 
carried out by the HF-H 2 S0 4 method. The formation of these insoluble pre- 
cipitates produced large uncertainty in the abundance determination because 
of the possibility of loss of trace elements by coprecipitation. 

The solution obtained above was then adjusted to 1.5-N HCl by 
repeated evaporation with concentrated HCl and finally by dilution. The 
evaporation with concentrated HCl must be repeated until removal of nitrate 
ions is complete. 

For the abundance determinations, appropriate isotopic spikes must 
be added at the beginning of the decomposition step in order to 'achieve best 
mixing of spikes with the samples. 

C. Separation of Silver -by Ion Exchange 

It has been proposed to separate silver from 0. 1-N HCl by an anion- 
exchange column. As mentioned above, it is difficult to obtain complete 
volatilization of lead and thallium. Thus, the procedure was slightly modi- 
fied to take this into consideration. Lead, thallium, and silver are ab- 
sorbed on a Dowex 1 X 10 anion- exchange column from 1.5-N HCl solution at 
this point. Preliminary experiments using radioactive tracer indicated 
that silver is quantitatively absorbed onto the column, and that it can be 
successfully eluted by 6-N HCl. Lead can be eluted by 0. 1-N HCl from the 
column before the elution of silver. Thallium is then eluted by dilute 
H 2 S0 4 saturated with S0 2 . 

Rare-earth elements, uranium, and thorium were observed to remain 
in the 1.5-N HCl solution in this step. The effects of interf erring ions 
such as Fe 1 "*”* in this step* have not been fully investigated. Thus, further 
experiments on this procedure will be carried out, and modification of this 
step may result. 

D. Separation of Uranium and Thorium by Ion Exchange 

Previously, it was planned to separate uranium and thorium from 
8-N HN0 S solution by an anion- exchange column. The direct conversion of 
1.5-N HCl solution into 8-N HN0 3 solution by addition of nitric acid and 
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heating resulted in the formation o.f insoluble residues, which might ab- 
sorb trace amounts of uranium and thorium. Thus, a series of experiments - 
will be carried out to find the most adequate way of conversion. Most 
likely, hydroxide precipitation and dissolution of the pr'ecipi'ta'te in 
several portions of 8-N HN0 3 will be- added to the procedure. 

Further studies of the isolation technique for uranium and thorium 
are -being carried out at present. The separation- of uranium from thorium 
and purification of each will be performed using procedures .such as the 
ones, described’- elsewhere in thia report (5,6).. 

E. Separation of Rare-Earth Elements 

Rare-earth elements are at this point in the 8-N HNO3 solution with 
only a few other elements. The rare-earth elements will be separated as 
a group by fluoride and hydroxide precipitations. The separation of rare- 
earth elements from this point is discussed in more detail elsewhere in 
this report (7). 


III. Summary 


The proposed schemes for the analyses of lunar .samples are summarized 
as block diagrams , i-n-<Fi^ures— 2~and— S-r 3 

The separation scheme for the mass^spectrometric elemental abundance 
determination is very similar to the scheme for the preparation of samples 
for the isotopic-composition determination of the elements in question,,, 
exce pt -tha-tr~the formen—inolud&s—the_addii^Lon~of^~iso.t i opi-e^s plkes-'a-t— appro — 
ypriaf- e-s-t-Ap s-. For the determination of lead and thallium abundances, spikes 
are added to an aliquot of the solution o'f the volatile deposit. 

Shie— eho^&* 5 ^^the_v.o.la.t.^^ chemical isolation 

of lead and thallium was^based on the fal-lotfxng advantages: (a) the smaller 

chance of contamination from^x^eagen-ts7 (b) simplicity of the separation 
scheme, and (c) merit of add-i-tf^ia-l^f lexibility, such as differential/ heating 
of the sample. The mal-n'^isadvantage&f^this volatilization method is that 
it is extremely^di-ffocult to attain complete"s'epaxation of lead and thallium 
from the^matr-rx7 and the yields in individual runs^re-^uiiknown. 

The investigation of the Pd 1Q7 -Ag- t -°~ 7 system in lunar samples has 
been discontinued tentatively because of extremely low concentrations of 
Ag in Apollo 11 samples , and the indications of a rather long 

time from the termination of r-process nucleosynthesis to the formation of 
cool meteoritic solids However, because s-process nucleosynthesis 

may have terminated later than / the Ir-procefss , and fractionations of non- 
volatile elements may have occurred earlier than the last cooling to the 
point of rare-gas retention, the s larch for extinct natural radioactivity 
of Pd 107 continues to appear worthphile. / 

The proposed scheme d-xscussed here is not yet completed and is 
undoubtedly^'^'ubjec.fe^to furt-^er/modif ications . Upon completion of various 
extentiuhs o^th^schexn^f tl^y will \be also applied to the analysis of 
me t edr i t e^sjtpvpT e s--as 1 


those returned from the moon in the Apollo 
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Diagram of the vacuum volatilization system. 
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Block diagram of chemical separation scheme for elemental 
abundance determinations 
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Figure 3. Block diagram of chemical separation scheme for isotopic-composition 
determinations. 
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MASS -SPECTROMETRY 
Lance P. Black and James M. Huey 


I. Additional Instrumentation 

« 

The Nuclide model SU-2.4 mass-spectrometer (90° deflection, 30 -cm radius), 
described in a previous report (l), has been upgraded during the- past year with; the 
addition of' several new systems. This section discusses the reasons for installing 
these systems j- a discussion of their electronic circuitry can be found elsewhere 
In this report (2')‘. ' ‘ - 

A. Digital recording system 

A digitalizing system, similar to that discussed' by Comps ton, Lovering, and 
Vernon (3 ), has been constructed by mating a Hewlett Packard model 2212A voltage -to- 
frequeney converter with a Monsanto model 1510A counter-timer and a Hewlett Packard 
model 5050B printing digital recorder. This allows an additional output mode besides 
the normal and expanded -scale circuits of the chart recorder. The system is 
particularly useful for the measurement of rapidly decaying signals, which cannot 
he satisfactorily accommodated by the expanded -scale circuit. As the digital- 
voltmeter does not utilize the chart recorder, it is free of any errors arising from 
non-linearity in this component. In addition to advantages arising from its high 
sensitivity and precision, the digitalizing system presents the data in a more 
convenient form than previously. 

B. Magnet -switching device 

•Until recently it was possible to alter the magnetic field during data 
collection only by continuously varying the magnet current. A control has now been 
installed which allows the magnet current to be alternately switched between two 
preset and adjustable values . More time is now spent on -peak rather than at the 
less important off-peak position; the latter is monitored at the beginning and end 
of each set of measurements' of- a particular isotope ratio. Data can new be collected ’ 
in less time and are less affected by periodic shifts in beam intensity than previously 
The combination of the magnet -switching device and the digital -recording system produce: 
data which are readily amenable to statistical analysis. 

C . Water pump 

A Tuthill l/4 h.p., 1800 r.p.m. pump is now used to recirculate cooling water 
to the diffusion pump. This has eliminated the need for the city water supply, which 
previously caused fouling of the safety cut-off switch after only a few days use. 

The system is now totally immune to water pressure failure in the mains supply,* power 
failures will cause the diffusion pump to shut off before overheating. 

D. Bakeout tapes 

Two insulated electrical heating tapes attached to a rheostate control have 
been wrapped around the source chamber for the purpose of periodically removing 
volatile -material from the- walls .of the- vacuum - system. Each tape measures 1X9° 
and has a maximum output of 288 watts. The bake -put system has improved the vacuum 
within the instrument and considerably reduced the time required for pumpdown after 
sample introduction. 
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E. Optical pyrometer 

A Pyro micro-optical pyrometer Model No. 95 (manufactured by the Pyrometer 
Instrument Co., N.J. ) has been fitted to a permanent stand opposite the filament - 
observation portal. The pyrometer is capable of measuring temperatures from 700 °C 
to 3200 °C on targets which may be smaller than .OOl" in diameter. With this 
instrument, filament temperatures can be monitored to approximately + 5 °C. This 
allows a precise control on the thermal condition of the source, a parameter which 
can affect measured isotopic ratios. Precise filament -temperature control is 
particularly critical in the case of lead (see Part II. A.) 

II . Techniques 

A large part of many scientific investigations is devoted to developing 
satisfactory analytical techniques. This is particularly true of mass -spectrometry. 
The analysis of meteoritic and lunar samples has, in many cases, demanded the 
analysis of smaller quantities of certain elements than are routinely measured. 

Thus, relatively new methods have been applied, some of which are still in the 
developmental stage. In particular, the techniques used for Pb and TT analysis are 
still being modified. 

A. Lead 

The low lead concentrations reported for preliminary analysis of Mare 
Tranquillitatis surface material (4) necessitated the development of a highly 
sensitive analytical technique. Following discussions with Drs. Shields and Barnes 

(5) of the National Bureau of Standards (Washington, D.C.), a method based on the 
work of Cameron, Smith and Waller (6) was adopted. We have since used this in the 
analysis of Apollo 11 and 12 lunar fines {J } and unpublished data), and is presently 
being used for meteoritic lead analyses (8). The technique involves the use of 
silica gel during the thermal ionisation. This matrix retards the evaporation of 
lead ions until higher temperature, where there exists a more favorable ionisation- 
to-evaporation ratio. 

The gel is prepared by a method similar to that of Cameron, Smith and Waller 

(6) . However, four heatings of gel with 5$> HNQ3 are used instead of two. Great; 
care is also taken during this stage to ensure that the gel does not dry out, as 
this leads to coarse-grained particles unsuitable for mass -spectrometry. The gel 
is stored under a shallow water layer (~ O.25") in glass jars with tightly -fitting 
plastic caps. 

After stirring and a one -minute delay, a 10 -pi aliquot of the overlying aqueous 
layer, containing gel particles in suspension, is taken and placed on a single 
filament (zone -refined, high-purity rhenium measuring .030" x .0012 " \ supplied by 
the Bembar Co. ) . The lead sample in % HNCb is added, followed by a small quantity 
(l pi) of H3PO4. Varying strengths of this acid are used, as some batches of gel 
emit more stably with weak acid (0.25 M), while others require stronger acid (2.5 M). 
The sample is dried by resistance heating using a current of about 1.8 A in air for 
approximately five minutes. 

\ 

■ After introduction to the mass -spectrometer ^ionization chamber, the filament 
temperature is raised to 800 °C in the space of about ten minutes. It is then 
adjusted upwards by about 70 °C every five minutes until 1170 to 1190 °C. At these 
temperatures there is negligible hydrocarbon contamination. Data are taken after a 
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stabilization period of 10 minutes. In seme cases satisfactory stability is 
not achieved after this time. When this occurs, the filament current is 
taken to approximately twice its normal reading for about a second, then returned 
immediately to its norma! value (suggested in a personal communication from 
S. Moorbath (9)). lata taken, on the NBS 981 and 983 standards suggest that this 
treatment does not alter isotopic composition (see Table I). 

Unlike the more conventional solid -source procedures employing sulfide (10.) 
or oxalate (11), the silica-gel method can be used for 'sample sizes in the vicinity 
of, and less than, one microgram. A one-microgram lead sample yields total ion 
currents of between 1CT 11 and icn 10 A. Ion signals of 10 -la A, which persist for 
more than an hour, can be obtained from less than 100 ng of Fb. Hence we have found 
it possible to determine the isotopic composition of our chemical blanks (~ TO ng). 

A further advantage is that the method produces no measurable isotopic 
fractionation (see Table l). This is particularly useful as it avoids the time- 
consuming and contamination -prone double-spike technique. The table shows that 
our average composition for UBS 981 is not significantly different from Catanzaro' s 
(12) value (at the 95/> confidence level). Our 206/204 and 206/207 ratios of 
UBS 983 also agree with those obtained from the KBS triple -filament work; our 206/208 
ratio, however, appears a little low. In the light of the combined results on both 
standards, this anomaly is difficult to explain by either contamination or isotopic 
fractionation. 

Ion beams are not as stable as those in the more common solid-source procedures 
but are generally quite satisfactory, especially as only one mass -spectrometer run 
per analysis is required. 

B . Neodymium 

The loading procedures for Ifd are relatively simple. The sample, in 5 ^ 

HHO 3 , is placed on two rhenium side -filaments and then converted to the oxide by 
resistance heating at about 1.8 A. Side-filament temperatures between 1300 and 1400 
°C are used for data collection; the center filament is set at 1450 to 1500 °C. 

Each isotope is measured relative to that at mass number 146. Total beam current for 
a 5-fig sample is about 6 x 10” 11 A. 

The data in Table II have been obtained from three complete mass-spectrometric 
analyses on commercial reagent neodymium* these are in reasonable agreement with the 
published compositions (13, 1.4). Currently analyses are being made to find the 
reason for the small discrepancies in some of the ratios. Particular attention is 
being paid to stricter time and temperature controls, as it has been found that the 
l42/l46 ratio, at least, may change by approximately 1$ during the early stages of a 
run. 


C . Samarium 


Loading procedures are the same as for Ifd, with the nitrate being converted 
to oxide on the filament. However, as with Nd, it is the singly-charged atomic ions 
that are used for analysis. A center -filament temperature of approximately I 65 O °C, 
used in conjunction with side -filament temperatures of around 1300 to 1400 °C, 
appears to give the best results. Each isotope is measured with respect to that at 
mass number 149. Ion beams are generally about twice as strong as those for Hd. The 
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reproducibility shown in Table III, which was obtained under 'a wide range of time 
and temperature controls, should be improved by experimentation under more rigid 
conditions . 

D. Thallium 

Continued experimentation with the triple -filament hydroxide technique 
reported by Huey (17) has shown it to be inadequate for the analysis of the low-ppb 
levels of the lunar materials. The single -filament silica -gel method used for 
lead analysis has been investigated,- but has also proved unsatisfactory. Although 
the technique does have the desired sensitivity, it. does not produce consistent 
isotopic ratios at the relatively low operating temperatures of less than 800 °C. 

The inability to precisely monitor these filament temperatures appears to lead to 
variable isotopic fractionation and/or variable amounts of hydrocarbon interference. 

A combination of triple -rhenium filament techniques with the silica -gel 
-method of sample preparation -provided more consistent results for 25-50 ng samples 
of thallium. However the low concentrations of thallium observed in many natural 
materials, including meteorites and lunar fines, makes a method for mass analysis 
of <5 ng samples desirable. 

Mass analysis of such small samples yielded consistently low 205/203 ratios, 

2-.386 > R > 2.220, relative to literature values of 2.590 (18, 21, 22, 25, 24) 
depending on filament current , time of heating, and' amounts of silica gel and 
phosphoric acid used. The presence of .a small peak at mass number 201 suggested that 
Re0+ ions (mass numbers 201 and 205) produced during heating with phosphoric acid 
may be responsible for a low 205/203 ratio. Accordingly, tungsten was substituted as 
a filament material. 

Tungsten -ribbon side filaments are loaded with silica gel, sample, and ~ 2 ul 
2.5-M;E3P04. Drying is accomplished by resistance heating (2.0-2.5A, 15 min.). 

During analysis the ionizing filament is heated to 1000° (~2A) and side -filament currents 
.are turned up slowly. The presence of hydrocarbons is observed at low temperatures, 
and data are collected after they burn off (side -filament T = 800°(~1.3A)). 

Table IV shows the results of several analyses obtained by this technique in 
comparison with some previous- results. There is an apparent increase in the precision 
of an individual analysis with this technique as compared to previous methods 
utilized in this 1 laboratory (18), In addition there is a rather large increase in the 
precision of replicate analyses. This can be attributed to the. greater sensitivity of 
the method (larger signal obtained) ^the stability of the signal, and the use of a 
digital readout system. 

This technique is being applied presently to the analysis of meteoritic thallium 
to determine thallium contents (using stable - isotope dilution) and isotopic 
composition (19) • 
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Table I 

Isotopic Analyses of N.B.S. Lead. Standards 


ms 981 

208/206 

206/20k 

207/206 

Analysis 68 

2.168k 

16.909 

. 9150 ; 9136 * 

102 

2.1697 

16.909 

.9169 

106 

2.1690; 2.1698 

16 . 95 k 

.9161 

127 

2.1627 

16.92k 

•9lkl 

168 

■ 2.16575 2.165k* 

16.916 

.9150 

175 

2.1632* 

16 . 916 * 

•9lk6* 

179 

2.166k* 

I6.9k8* 

• 91 k 9 * 

Average composition 

2.1665 + .ooio§ 

16.925 + *°° 7 § 

.9150 + .000k§ 

Ef.B.S. Value ( 12 ) 

2.1681 

16.938 

.9lk6 ' 

NBS 985 

206/208 

206/20k 

206/207 

Analysis 115 

72.95 

2622 

lk.oko 

185 

73.25 

2835 

lk. 038 

187 

73 . 15 - 

27985 2595 * 

Ik.Oki 

190 

73.22 

2791 ; 28 e 9 s 2755 * 

lk .033 

Average composition 

73 . 1 k +_ , 07 § 

27k8 + k7§ 

' lk,038 + .002§ 

N.B.S. Value ( 12 ) 

73 .k 

2698 

1 

lk.ok 

*Data sets taken at beginning and end of run, 
★Data taken at three stages during run. 

respectively. 



*Data taken after sample was briefly subjected to high filament temperature. 
§Uheertainties quoted' represent the standard deviations of the means . 


Table II 


Isotopic Analyses of Commercial Reagent Neodymium 



1kg 

lk 3 

ikk 

ik 5 

Xk6 

lk8 

150 

Analysis 92 

27.08 

12.16 

23.77 

8.30 

’ 17.25 

5.78 

5*66 

97 

27.09 

12.20 

23.80 

8.30 

17.20 

5.76 

' 5.65 

10k 

27.05 

12.18 

23.77 

8.30 

17.23 

5.79 

5.69 

Inghram et al.( 13 ) 

27.13 

12.20 

23.87 

8.30 

17.18 

5.72 

5.60 

White et al.^lk) 

27.3 

12.32 

23.8 

8.29 ' 

17.10 

5.67 

5.56 
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Table IU 

Isotopic Analyses of Commercial Reagent Samarium 



144 

147 

148 

M 

150 

152 

154 

Analysis 20 

3.09 

15.00 ' 

11.24 

13.81 

7.38 

26.75 ■ 

22.72 

96 

3.01 

15.00 

11.26 

13.86 

7.43 

26.76 

22.68 

100 

'3.11 

15.07 

11.30 

15.86 

7.38 

26.67 

22.61 

101 

3.08 

15.02 

11.26 

13.84 

7.39 

26.70 

22.71 

Inghram et al.(l5) 

3.16 

15.07 

11.27 

13.84 

7.47 

26.65 

22.53 

Collins' et al.(l 6 ) 

3.15 

15.09 

11.35 

13.96 

7.47 

26.55 

22.43 


Table IV 

Isotopic Analyses of Reagent Thallium 



Sample Size 

'205/203 

Analysis 248 

< 5 ng 

2.389 + .004 

252 

< 5 ng 

2.388 + .003 

253 

< 5 ng 

2.382 + .002 

257 

< 5 ng 

2.386 + .001 

259 

< 5 ng 

2.390 + .003 

Weighted Avg. 

< 5 ng 

2.388 + ,002 

Avg. of Ostie's Replicate Anal. 
Normal Tl* (20) 

Normal Tl (21, 22, 23, 24) 

( 3 . 8 ) 40-80 ng 

2.389 + .005 
2.36 + .02 

2.390 


★ 


correction for electron. multiplier applied 
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II. D. 4. 

CHEMICAL PROCEDURES FOR PURIFICATION OF URANIUM AND 
THE ALPHA-SPECTROSCOPY OF URANIUM 

Nancy -Kan 


I . Introduction 


Tektites are glassy objects found scattered on the Earth's surface 
in discrete regions up to thousands of kilometers apart. Their origin is 
unknown. There are two major hypotheses, lunar origin and terrestrial origin. 
Schwarcz (1) has demonstrated similarities in the chemical composition of 
tektites and soils, and assumed that tektites are formed on the earth from 23g 
soils by the impact of some extra-terrestrial object. In this case, U /U 
isotope ratios differing from unity in tektites might provide strong con- 
firming evidence of this hypothesis (2,3). The U 2 4 /U 238 isotope ratio will 
be an important clue to the origin of tektites. 

If the young tektites showed a U 234 /U 238 ratio below unity and old 
tektites did not, this would prove that tektites are terrestrial in 3 ori|in. 

So we have undertaken to determine as accurately as possible the U /U 
ratios in tektites.. This will require higher precision than in previous work, 
~0.1%. This report describes a chemical procedure of isolation and purifica- 
tion of traces of uranium to determine the U 234 /U 238 activity ratio of alpha- 
particle spectrometry. 

II. Chemical Procedure 


The techniques of decomposition of silicate rocks and chemical separa- 
tion are those developed by T. P. Sarma (4). Vassiliki Tsomi (2) followed the 
same procedure. Thomas J. Thomas (3) modified it by using a thenoyl-tri- 
fluoro-acetone (TTA) extraction to deposit uranium onto a planchet, which 
has also been used by Ku (5) . The TTA method was found (3) to yield samples 
with better alpha-spectrum resolution than samples prepared by electrodeposi- 
tion. This may be due to the more rapid deposition in the TTA method with 
correspondingly less granular crystal structure. Also, Thomas changed the 
planchet material from platinum to titanium ( 6 ) . The chemical procedure Is 
outlined in Figure 1. 

A. Decomposition of Silicates 

After treating the 10-g sample with HF, H2SO4, and perchloric acid 
and leaching with HCl as in the procedure of Sarma (4),' we still have some 
residue which cannot be dissolved. For dissolving the residue, instead of 
using a sodium fusion, to avoid contamination, a mixture of equal volumes of 
HF, HCl, and HNO3 is used in a teflon beaker and boiled with the residue. 
Several portions of HCl are added and evaporated to get rid of the HF and 
HNO 3 . The solution is heated under a heat lamp for several hours to reduce 
the volume to 500 ml. The solution is cooled and 5 g of sodium bromate is 
added in small portions in order to oxidize U(IV) to U(VI). Then the solution 
is neutralized with concentrated ammonium hydroxide until it is slightly basic . 



The solution is heated under a heat lamp for a couple of hours . The hydroxide 
precipitate is centrifuged and washed twice with water and then evaporated to 
dryness . 


B. Ion-exchange Separation 

A polyethylene anion- exchange column containing 30 ml of Dowex 
1-X4 is set up and washed first with 150 ml 1-M HC1 and then with 150 ml 
10-M HC1. The dried sample obtained in the co-precipitation step is dis- 
solved in 600 ml 10-M HC1 and transferred to the ion-exchange column. The 
column is washed wi-th 4-column volumes of 10-M HC1, when all the major 
elements and radioelements except Fe(III), Pa(V), U(VT), and Po(IV) are 
eluted from the column. Leaching the column is continued with a solution 
of 150 ml 9-M HC1 plus 1-M HF to remove Pa(V) . Finally uranium is stripped . 
off from the column with 300 ml lvM HC1. The .flow rate of the 'elution is 
about 8-10 drops/min. 


G. Di-isopropyl Ketone Extraction 

The uranium fraction containing Fe(III) and Po(IV) is evaporated 
to dryness and taken to fumes with concentrated HNO3 to eliminate the last 
traces of fluoride ion.. The residue is dissolved in 25 ml of 6-M HC1. The 
solution is extracted with di-isopropyl ketone by the countercurrent extrac- 
tion method, using three separatory funnels, when Fe(III) and Po(IV) go • into 
the organic phase. The aqueous phase containing uranium is evaporated to 
dryness with 50 ml concentrated HNO3 to decompose any organic matter. 

D. TTA Extraction 


The thenoly-trifluoro-acetone (TTA) method has been used' by Ju (5) 
to deposit uranium onto a planchet. The dry residue is taken up with 0.2. ml of 
0.1-M HNO3. The purified uranium is stripped into the TTA phase in a 0.4- 
M. TTA-benzene extraction after the pH of the aqueous phase has been adjusted 
to 3.0 by using a pH 3 buffer. The extractions are done by countercurrent 
extraction with equal volumes’ (1 ml each) of TTA and uranium solution. 

-E. Preparation of Sample for Alpha Counting 

The benzene-TTA-U phase is evaporated onto a Pt or Ti disk which 
is washed first with acetone, then with concentrated HNO3, and finally with 
distilled water. The disk is then flamed to burn off the organic matter, 
giving a thin source of uranium oxide. Another ‘method of preparation of 
samples is that of direct evaporation on a Ti disk of a pure uranium solution 
with the aid of a spreading agent — * tetraethylene glycol (TEG) (7) . About 
1 or 2 drops of TEG are dropped on a Ti disk, and the uranium solution is 
then transferred to the disk with a capillary dropper. The solution is 
evaporated slowly under the heat lamp to dryness. The disk is then flamed 
and the uranium is left as a uniform thin. film. 
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HI* Alpha— Particle Spectrometry 

The alpha spectra are obtained by an alpha system consisting of 
a 450-mm ORTEC detector. Instrument and Development Products preamplifier; 
amplifier, and a Nuclear Data 256-channel analyzer ( 8 ). 


A typical spectrum of uranium from uranium nitrate ( 1102 (^ 3)2 
6 H 2 O, Baker analyzed reagent grade) is given in Figure 2, This sample 
contains U 34 » D 38 , U 23 , and also t-races of U 236 . The U 236 is present 
as a result of intense neutron irradiation, evidently in a nuclear reactor 
from whose spent fuel elements this uranium was recovered. The facts that 
the U /U 23 ratio is not unity and the pr 23 s /u 238 ratio is below that in 
natural uranium also indicate that the -sample- is not of natural uranium. 
Because of the deficiency of U 234 , the sample can be used for isotope dilu- 
tion analysis. 


The U 234 /U 238 , U 235 /U 238 , U 236 /U 238 ratios are calculated by taking 
ttlS i ? 38 8ral ° f the al Pk a spectrum. (Figure 2) from channel 96 to channel 132 
’ c ^ ianne ^ -f-33 to channel 155 for U 235 , channel 156 to channel 165 for 
U 23 s’ and channel 166 to channel 2D2 for U 234 . The widths of the -U 2 3 4 and 
U peaks are about 5 channels (-40 keV) FWHM. A comparison of the curve 
in Figure 2 with the synthetic natural uranium spectra of Thomas (3) , shows 
that the resolution is somewhat better than the 40-keV FWHM curve and de- 
finitely worse than the 20— keV FWHM. curve. The instrumental resolution is 
thus about 35 keV. 


IV. Quantitative Test of Procedures 

The counting yield at different shelf positions can be calculated 
by using a uranium’ nitrate solution and assuming that the specific activity 
of U 238 is 734 disintegrations/mg/min (Table 1). 

The data have been collected by using a uranium nitrate solution 
to check the chemical yield and the reproducibility of the whole procedure. 
This is accomplished by performing only the last step of the procedure in 
the first run, and going backwards and adding one more step, proceeding to 
the beginning, for. each successive run. For, example, first only the TTA 
extraction is done, next the di-isopropyl ketone extraction followed by the 
TTA extraction, and so on. The data from successive addition of steps shows 
very good agreement in the IJ 234 /u 236 ratio (Table 2). The chemical yield is 
based on a sample of uranium solution directly deposited on a Ti disk by the 
TEG spreading method (Table 2). 


V. Plans 


Because of the low-energy tailing of alpha spectra, low instrumental 
resolution from thick samples, and the necessity of getting higher precision, 
it is necessary to use a computer analysis of data. This was begun by 
Thomas (3) and is being developed further by Trivedi (9) . 
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By using a U tracer through the isotope dilution method, the 
chemical yields at very low uranium concentration can be tested. The goal 
of this whole procedure to try to get the uranium content of tektites by 
isotope dilution and the U 234 /U 238 activity ratio in tektites will thus be 
fulfilled. 
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Table I 

COUNTING YIELDS OF 250 yg URANIUM DEPOSITED ON* Ti 
WITH TEG (SAMPLE U-15) AT DIFFERENT POSITIONS OF THE DETECTOR 



Activity Ratio 

Counting 
Rate ‘for 
U 238 . Peak 

Counting 

Yield 

Position 




Third' Shelf • 
(0.937 inch) 

0.232 

0,0178 

0.0077 

709/hr 

6.4% 

Second Shelf 
(0.687 inch) 

0.241 

m 

-0.0093 

1109/hr 

10.1% 

First Shelf 
(0.437 inch) 

m 

0.0198 

0.0090 

1872/hr 

17.0% 






















Table -II 


CHEMICAL YIELDS FOR PROCEDURE STARTING 
AT VARIOUS STEPS 


Sample 
(500 yg 
Uranium)- 

Procedure 

Activity Ratio 

Counting 
Rate for 
U 238 Peak 
' (3rd Shelf) 

Chemical 

Yield 

u 234 /u 238 

3 5 j2 3 8 

|j2 3 6 3 8 

U“4 

Starting with 
TTA extrac- 
tion 

0,239 

0.0189 

0.0082 

850/hr 

60% 1 

U-5 

Starting with 
TTA extrac- 
tion 

0:236 

6.0189 

0.0079 

808/hr 

58% 

U-22 

Starting with 
TTA extrac- • 
tion 

0^239 

0.0243- 

0.0081 

1045/hr 

' 74% 

' U-r-20- 

Starting with * 
di-isopropyl 
ketone ex- 
traction 

0.233 

0.0193 

0.0081 

579/hr 

' 

41% 

, U-14 

Starting with 
ion-exchange 

0.233 ■ 

• 0.0184 

0 : .0074 

- 5 33 /hr 

- 38% 

U-16 

Starting with 
ion-exchange 

0.244 

0.0171 

0.0067 

• 4 25 /hr- 

30% 

U-17 

Starting with 
hydroxide 
precipitation 

0.232 

0.0257 

0.0118 

394/hr 

28% 



Figure 1. Chemical Procedure for Uranium 
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Figure 2, Alpha Spectrum of Commercial Uranium, 
Sample U-16 (See Table 2 for prepara-' 
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II. D. 5 



alpha-isotope-dilution analysis of uranium and thorium 


10 

* -X 


Brij M. P. Trivedi 


The present work aims at quantitative determination of uranium 
and thorium in natural samples, with special interest in meteorites and 
lunar materials by the radiometric alpha-isotope-dilution technique using 
U 232 and Th 228 as spikes. 


I. Chemical Procedure 


A. U 232 and Th 229 Tracer Solution 


U 232 solution was purchased from General Nuclear, Inc. A 
standard solution of 0.01 pc/ml of U 232 was prepared by T. J. Thomas from 
the original standard. A portion of this solution was further diluted to 
give 90 dpm/ml of U 232 . About 0.1 ml of this solution was used as spike 
in the development of the chemical procedure. The activity ratio of 
(U 232 /Th 22S ) is about 6. 

B. Dissolution of Silicate Samples 

A weighed quantity 2-5 g) of the silicate sample is taken in 
a teflon beaker and spiked with 0.1 ml of U 232 standard solution. The 
same amount of U 232 solution is evaporated on a platinum planchet and 
counted separately to check the chemical yield. About 25 ml of hydro- 
fluoric acid and 25 ml of nitric acid are added to the sample. The 
solution is evaporated to dryness and the HF + HN0 3 process repeated 
three times. The residue is dissolved in concentrated HC1. In .all cases 
a clear solution is obtained. 

C. Group Separation 

Chemical procedure for separation and purification of U and Th 
has been developed by many investigators (1, 2 , 3 ,4, 5, 6, 7, 8). Briefly, 
the solution obtained after dissolution of silicate material is evapo- 
rated to dryness and taken into 1.8-MHC1. The solution is passed through 
a Dowex 1 X 10 anion-exchange column (200-400 mesh, Cl form) 20 ml in 
volume. The column is washed with 1.8-MHC1 (6 column volumes) when U, 

Th and rare earths are separated from the column and silver sticks to the 
column. The eluate is dried and taken into8’-MHN0 3 and passed through 
Dowex 1 x 10 column (20 ml volume) . Rare earths are separated from the 
column by washing it with8-MHN0 3 . Thorium is separated from the column 
by washing it withlO-MHCl and finally uranium is separated by washing 
it with 0.1 -M HC1. 
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Jr 'C VI 

gjy f PQTLf-icatibn o3: Thorium 


The thorium fraction is evaporated to dryness and taken into 
8-MHCl. The solution is passed through Dowex 50 WX4 cation-exchange 
column (ionic form fl+, 50-100 mesh, '20 ml volume) . Thorium sticks to 
the column. Thorium is desorbed from the column by washing it with 
0. 5-Moxalic acid. Oxalic acid in the solution containing' thorium is 
decomposed by heating it with concentrated nitric acid. The solution 
is evaporated to dryness and taken into 10-M HCland passed through 
Dowex 1x4- anion-exchange (Cl" form, 100-200 mesh) column (20 ml 
volume). Washing the column with four ( column volumes of 10-M HC1 removes all 
thorium. Pinal purification of thorium is done by extracting it with 
0.4-M solution of TTA (thenoyltrif luoroacetone) in benzene. Thorium 
solution is evaporated to dryness and taken into .approximately half ml 
of 0. 1-MHNO3. Equal volume of 0.4-MTTA solution is added and equili- 
brated with the aqueous solution. Thorium comes out in the organic 
phase and is deposited by evaporation on a platinum planchet. The 
platinum planchet is.. put to flame, when a very thin source of thorium 
is obtained. ** — " 


E. Purification of Uranium 


The uranium fraction is evaporated to dryness and taken into 
6-M HC1 (~ 15 ml) . The solution is equilibrated with an equal volume of 
diisopropyl keton'e. Uranium remains in, the aqueous phase and is free of 
iron. The uranium fraction is heated with concentrated HN0 3 to decompose 
any. organic matter and evaporated to dryness. The residue is taken into 
10-MHCl and passed through Dowex 1x4 anion-exchange (Cl” form, 100-200 
mesh) column (20 ml). Uranium sticks to the column and is desorbed from 
the column by washing it with 0.1-MHCl. The solution is evaporated to 
dryness- and taken in about 0.2 ml of 0.1-MHN0 3 . A buffer solution 
(prepared from buffer tablets, Coleman) of pH 3 (~ 0.5- ml) is added to. 
make the pH of the solution about 3. The solution is now equilibrated 
with equal volume of TTA (04-M in benzene), when uranium comes out in 
the organic phase. The organic phase is withdrawn and evaporated on a 
platinum planchet to deposit the uranium uniformly. The planchet is 
heated -red hot, giving a thin source of uranium. 


II. Alpha-Particle Spectrometry 

The samples are counted in a surface barrier detector (Model No. 

E -022-450-300) , with active area 450 mm 2 and sensitive depth 300 microns, 
purchased from 0RTEC. The resolution of the detector is 35 keV and the 
counting yield 18%. 


III. Preliminary Experiments 

Preliminary experiments have been carried out with tracer alone 
(U 232 anc j Th 22S solution) to check the procedure. Figures ( 1 ) and ( 2 ) 
show the' alpha spectra of U and Th, respectively, separated from a tracer 
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solution. Further experiments have been carried out with pure silica 
glass with spike added (blank run) and rock sample' (Granodiorite ; GSP-1) 
obtained from the U. S. Geological Survey. The chemical yield for each 
case was found to be about 50%. Figure (3) shows the alpha spectrum of 
uranium separated from the rock sample. The thorium sample is being 
purified for counting-. 
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Figure 1. Alpha spectrum of U 232 separated from a tracer 

solution containing U 232 and Th 223 and its decay 
produc ts . 
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Figure 2. Alpha spectrum of Th 228 separated from 
a tracer solution containing U 232 and 
Th 228 and its decay products. 
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II. D. 6 



LEAST-SQUARES ANALYSIS OF EXPERIMENTAL' ALPHA SPECTRA' 

I 

I 

Brij M. P. Trivedi 


For precise determination of the U 234 /U 238 ratio of uranium from 
natural materials (tektites, rocks, lunar samples, etc.) and for alpha- 
isotope-dilution analysis of uranium and thorium using U 232 or U 235 and_, 
Th 228 or Th 23 °, respectively as spikes, a rigorous least-squares program 
has been developed to analyze the alpha spectra of these elements. 

Observed spectra of the samples show considerable low-energy 
tailing because of energy loss from the finite thickness of the sample, 
back scattering of alpha particles from the supporting material, and 
the range of angles at which the particles pass through the window of 
the semiconductor detector in high counting geometry. The- alpha spectra 
are therefore represented by a modified Gaussian function of the form: 


I = I exp 
o 


[ (E •- E q )/W] M i ' 


L 1 + [1 - sign (E‘ - B ) ] [S (E - E ) i** 2 - 


[1] - 


where 

I = intensity in each 'channel due to an alpha group 

I -= peak intensity of an alpha group 

E = Energy 

E q = peak energy of an alpha group 

W = resolution parameter 

sign (E - E ) = 1 if E > E 

o o 

= -1 if E < E 

o 

- a if E = E 

o 

S = skewing parameter. S is derived from another parameter 
Q," called' tailing parameter, equal to ratio of peak 
height to low energy asymptote of tail. 


Q 


exp. (l/2W 2 S 2 ) 
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• *\ 


Equation [1] was first used by Huey at al. (ll, keeping - M 2 = 
2. It was modified by Thomas (2) by keeping. M L = 2 but allowing M 2 to be 
variable. Equation [1] is a more generalized form of these equations. 


Because the spectra are taken in a multichannel analyzer, the 
energy E is expressed in terms of the channel number of the analyzer: 


where 


and 


E = C x + 

= zero channel energy intercept 

(bias) of the multichannel analyzer 

C 2 = energy increment per channel 

X = channel number 


The parameters to be adjusted in equation [1] are the intensities 
of the individual alpha emitters (I 232 , I 2 34 > I 2 35 , ^ 236 > ^ 238 ’ or ^- 22 a> -^230 
l£ 32 ) for uranium or thorium samples, respectively, and the spectrum 
parameters W, S, C lt C 2 , M and M 2 . Henceforth, the intensity and the 
spectrum parameters will be collectively referred to as parameters. The 
peak intensity of each group is derived from its known branching fraction 
and the total intensity of the emitter to which- it belongs. 

For the least-squares analysis, the values of the above-mentioned 
parameters are first approximately guessed. The function I (I 232 » I , 

1 235 j * 236 > ^ 23 S> W, S, C l9 ^i» H 2 ) is -expanded about these parameters 
in a Taylor series, keeping only the first terms: 

1 - (i) ° + + + • • • • t21 

(I) is the value of the function I for the initial (guessed) set 
° ‘ / ai \ 

of parameters. - The derivatives { ^ ) , etc. are evaluated numerically 

\ d - L 23^/° 

using the following equation and its analogs for other arguments: 




-*-^232 -*■ ^-232> l-234> * * * ’) ~ 1(1 


232> J -234 


*) 


AI 


[3] 


232 


AI 232 is 3 sma ll increment (for the present case it was 0.1 % of 
the original value) given to the parameter I 232 .. 

The condition of lease-square fitting adopted is (3): 


2 = E 


i -v 


I. - (I . 
x obs 


). 

1 


min. 


[4] 
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Summation is. over all the observations, i. a £ 2 = (Iobs)i (see 
for example Bevington 1969). For the present calculations <?i 2 ;Was set . 
equal to (Ical)i- This change was done to avoid the division .of the 
function [Ii**(Iobs) iJ by zero when (Tobs)l equals zero-. This change is- 
justified for the cases where Ij; is not much different from (l 0 bs^i - 
However., when the difference [Iir (I 0 b s ) j.] is large, o ^ 2 = C-Iobs^i should 
be used and a± 2 should be put equal to unity for the cases where .(lobs) 3 . 
equals zero. Other weighting factors have also been suggested in the 
literature (4,5,6). 

Equation [4] can be written as: 



x?here e. = [(I) -I,] for the channel i.- 

1 1 o obs 

Differentiating equation [5] with respect to AI 232 we get: 



•Similar equations are obtained for other parameters. These 
equations are linear in parameters (AI 232 , AI 23 4 ., etc.) and can be solved 
for them by psing matrix algebra. 


The goodness of fit at each step is judged from x^ s value defined 


s _ X* - _L 


(I - I ,.J ; 


N-n N-n j [ o '. 2 obs" 


where 


and 


N = total number of observations 
n = total number of parameters 
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The computer program has been written in ALGOL for the CMU 
computer center UNIVAC 1108 computer.' -It ac'cepts as input' the approxi- 
mate values of the initial phramters ^2345 ^236> I233, SjG^., 

C a , M a ) and the observed ' intensity in each channel.' It calculates 
total intensity in each channel i compares with the observed 'values finds 
an improved set of parameters , ' and goes' through several cycles to deter- 
mine the best 'set; So far each program was run only for five cycles, 
and the parameters were taken' for lowest value of x 2 . The program gives 
as output the calculated intensity in each channel, v the improved set of 
parameters and their standard derivatives, and the value of Xy 2 for each 
cycle of repetition. The method seems to converge in about 5-6 cycles if 
the initial parameters are properly guessed. 


Results 


The program has been tested by using experimental 'data on the 
measurement of uranium from tektites . The samples for these data were 
prepared by Tsomi (7) and counted by Thomas (2) . Two such fitted spectra 
are shown in Figures ‘1 and 2. For Figure 1 the observed spectrum was 
taken for a single run. For Figure 2, the same sample was counted at 
different times and the activity in different channels were added to give 
one spectrum. This might result in a poor spectrum when there is a large 
drift. Figure 1 shows the alpha spectrum of U from a moldavite [sample 
T-5-I-b; Tsomi (7), Thomas (2)J and Figure 2 shows the alpha spectrum of 
U from a philippinite [sample T-24-V,’ Tsomi (7)-, Thomas (2)]. Table I 
shows the initial set of parameters and the successively adjusted para- 
meters for each 'cycle for Figure 1. Also shown are the values 'rif X 2 for 
each cycle at the bottom. 

In Table II are given the derived .activity ratios of (U a34 /U a3S ) 
and (U a35 /U 238 ) for the two tektites. For comparison these ratios are 
given for equilibrium natural uranium. The isotopic ratios in tektites 
seem close to natural uranium isotopic ratio considering their relatively 
large uncertainties. More samples are being analyzed in this laboratory 
for isotopic studies (8). The method will be applied to isotope-dilution 
analysis of uranium and thorium from silicate materials (special interest 
is in lunar samples) .(9,). 


Discus sion . 


For a good fit, the value of X 2 should be close to unity (a 
value between 1 and 1.5 is reasonable) V For the "sample T-24-V, x v 2 
seems to be significantly different from unity. This might be due to 
the following reasons: (1) The method of analytical expansion of the 

function. I converges very fast only near the point of minimum of X 2 atl 4 
not from far off points. In other words, if the initial set of para- 
meters are not close to the true parameters, the method cannot be. expected 
to converge with any reliability. This difficulty can be overcome by 
using another method, which converges fast when. away from x 2 minimum but 
very slowly near x 2 minimum, to derive the initial set of parameters (3). 
In this method the function X 2 is minimized with respect to each parameter 
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separately keeping all others constant. The parameters giving the 
minimum value of y 2 can, then be used as the initial parameters in ^ 
the analytical expansion method. (2) Use of weighting factor (Wi = 

1 f l 1 

= - — ) may not be the best. Other weighting factors ( Wi = -vz vg; 

leal . \ UobS'' 

— r jwill be tried in .the future.- (3)- Incorrect nuclear data 

(I obs X 

(especially branching ratio of U 234 - alphas-) . This discrepancy is obvious 
from both Figure 1 and 2. (-4). .Variable resolution (functions of energy)-. 

The computer program developed is still in the initial stage. All 
the correction factors will be introduced, so that it can be used with full 
reliability. 
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TABLE I 

, ’ * , r 

ADJUSTED PARAMETERS FOR SAMPLE NO. T-5-I-b 


Parameter 

Initial 

Input 

Estimate 

- 

Adiusted Value After Cycle 

No. 

Standard 
Deviation 
± • 

1 

. 2 

3 

4 

5, 

W 

0.024 , 

0.021 

0.020 

0.020 

0.021 

0.021 

0.0000 

S 

12.30 

15.38 

18.81 

' 19.89 

19.64 

19.55 

1.61 

Ci 

3.360 

3.357 

3.355 

3.354 

3.353 

3.352 

0.0000 

C 2 

0.008 

0.008 

0.008 

0.008 

0.008 

0.008 

0.0000 

M x 

2.0 

1.92 

2.09 

2.22 . 

2.25 

2.26 - 

0.18 

m 2 

. 2.0 

•1.82 

' 1.90 

2.06 

2.11 

2.12 

O'. 08 

1 234 

260 

243.4 

250.4 

250.7 

245.6 ' 

245.5 

11.7 

235 

14- . 

14.34 

14.90 

14.63 

14.17 

14.07' 

2.13 

1 238 

260 

218.7 

223 .0 

223.2 

218.9 

218.8 

10.4 

ul 

2.63 

1.54 



2.97 

1.72 

1.46 

1,36 

■ • 


TABLE II 

ACTIVITY RATIO OF URANIUM ISOTOPES IN TEKTITES 


Sample Number 


No. of Iterations 
• W 
S 
C, 


l 234 / -*-238 




5 



4 


0.021 

± 

0.0000 

0.022 

+ 

0.0000 

19.55 


1.61 

14.97 


0.32 

3.352 

+ 

0.0000 

3.326 

± 

0.0000 

0.008 

+ 

0.0000 

0.008 

± 

0.0000 

2.26 

+ 

0.18 

1.85 


0.039 

2.12 

+ 

0.08 

1.72 

+ 

0.023 


245.5 


+ 11.7 


14.07 + 2.13 

218.8 + 10.4 


1.122 + 0.075 
0.0643 + 0.0099 


2012.0 

86.9 

2015.1 


+ 32.03 

± 5 - 64 
+ 32.01 
5.94 


0.997 + 0.022 

0.0458 + 0.00302 
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Figure 1. Experimental (points) and fitted (curve) spectra of 
uranium from a moldavite sample (T-S-T-h) 
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-II. D. 7. 


DESIGN AND 'CONSTRUCTION OF A 'POSITION-SENSITIVE PROPORTIONAL COUNTER 


John R. Yocum, Jr. and John F. Monahan 


I. Introduction 


The objective o£ this project is to design and build a propor- 
tional counter which will record the location of an event in the counter 
as well as its magnitude. In addition, we hope to make the counter light 
enough to be used in unmanned spacecraft. 

• i * ' • . 

The basis for our interest in position sensitive proportional 
detectors is that they would enormously decrease the complexity and in- 
strumentation of any proportional array and would tend to increase its 
sensitivity.. In a normal array, the smallest increment- of position is 
the size of the individual detectors. However, the. position increment 
of' a position sensitive device is dependent -only on the uniformity of 
the detector, a factor which can easily be controlled. There is a prac- 
. tical minimum size for proportional detectors at this time; but, theoret- 
ically at least, there is no limit to the accuracy of a position sensitive 
device. 

Another valuable point is that in order to achieve the accuracy 
now obtainable, many small detector.s, each with its own voltage supply, 
preamplifier, amplifier, and recording equipment, .must be used. If a 
large-=position sensitive device were used in place of each line of small 
detectors, the same or greater accuracy in positioning and recording the 
pulses could be obtained, at a much lower cost in. instrumentation,, weight, 
and power. ’ . . 

Two types of position-sensitive .proportional detectors have been 
considered. 

The first, a semiconductive detector, has been in use. for several 
years; however,, it has many serious limitations. It has very good accuracy 
on higher energy pulses, but it .also has an extremely high noise factor 
which blanks but lower' energy pulses. The size of one of these detectors 
is not very large, so an array of them would require the same amount of 
instrumentation as any other type of detector. The use of semiconductive 
detectors on unmanned spacecraft would be very limited, because detectors 
of this type require a low temperature to operate and are very sensitive 
to temperature change. Due to' the -latter consideration, cryogenic equip- 
ment must be used, adding to weight problems. 

-The other alternative,- a resistive proportional tube, is the main' 
focus of our efforts- at this time. In this type ’of' proportional* tube, 
the conductive center wire is to be replaced by a wire with a resistance 
of approximately twice the input impedance of one of the preamplifiers 
used with it. 
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The tube is to be operated by placing a high positive potential 
on the center wire from both ends of the wire. An amplifier and pream- 
plifier would also be attached to each end of the wire so that the circuit 
would be similar to Fig. 1. 
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Figure 1 


Any pulses occurring nearer one end of the tube -than the other, 
would be larger on that end because of the lower resistance path. By 
recording the sum of the two pulses, the magnitude of the event could be 
determined. The ratio of the two would indicate where the event occurred. 

The resistive proportional tube would have a wide range of appli- 
cations. Presently our interest extends mainly to cosmic X-ray mapping 
from unmanned spacecraft. Other uses in fields such as Medicine could 
also use the resistive proportional tube for cobalt therapy and organ scans. 


II. Experiments 

A. Resistive Proportional Tube Simulation 

In the resistive proportional tube simulation, a normal propor- 
tional tube and preamplifier were attached to the center tap of a potentio- 
meter and two amplifiers were connected to the other leads, as in Fig. 2. 
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The proportional tube used for the simulation was filled with a 
90% Kr - lo% CH 4 mixture- at 73 .cm pressure. Its operating voltage was 
1800 volts DC. Since the tube was especially designed for x-rays, we 
used a Co 57 source. 

By varying the setting of the potentiometer, the amount of signal 
reaching each amplifier could be altered. If the total resistance of the 
potentiometer is twice the input impedance of each amplifier and if the 
potential is assumed' to be distributed according to resistance only, the 
signal amplitude should vary from the normal value at one end of the 
scale to 4 value at the other. In the middle, each amplifier should re- 
ceive the normal pulse amplitude. - 

Two display techniques have, been developed. One places each am- 
plifier output on a beam of a twin-beam oscilloscope. Differences in 
pulse height and shape are very clearly indicated. The other method places 
one output on the horizontal deflection plates of an oscilloscope and the 
other output on the vertical. In theory, the result should be a sloping 
line whose length corresponds to the size of the pulse and whose slope 
corresponds to the ratio between the' outputs. However, due to inequities 
in the two amplification systems, noise and similar factors, the actual 
picture is very complex. Reduction of noise and more careful adjustment 
of the amplifiers' should alleviate this difficulty. We expect both of 
these display methods to evolve in ensuing research. 

• B. Production of High-Resistance Center Wire 

The high-resistance center wire necessary for the position-^sensitive 
porportional tube must have a -reasonable length -and strength and must, have 
a low- temperature coefficient of resistance. 

Two techniques of making" the high-resistance center wire are being 
considered: carbon films and vacuum deposition of metal films. 
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The carbon film is the first .to be considered and is the more- 
promising of the two methods. Carbon films are not as expensive to 
make, and they are more easily applied than metallic films. The thick- 
ness - hence, the resistance - of the film is easily adjusted, and it 
will not react with proportional counter gases such as argon-methane. ' 

Carbon films also have drawbacks, however. The primary difficulty 
is their low adherence to a glass or quartz fiber. Since glass or quartz 
will probably be used as the substrate for deposition,, a method of making 
the carbon film adhere to the substrate must be found if it is to be used. 

The actual deposition would probably be done from a gas; however, 
the details of this operation have not yet been worked out (1,2). 

If the carbon film proves unworkable, a nickel -chromium alloy 
deposited on a glass or quartz strand will probably be the second altern- 
ative chosen. The film would be vacuum deposited on an Si0 2 strand 15 cm 
in length to a depth such that the resistance would be approximately 
2 X 10 s & with a temperature coefficient of resistance of < . VL. The 
ratio of nickel and chromium would probably be 4:1 (3-6). 

Two methods of producing the 'film appear suitable. 

Flash evaporation calls for vaporizing small amounts of' the metallic 
powders in the correct ratio in very short periods of time. The best way 
would be to pour a slow stream of the powder onto a tantalum strip heated 
to about 1600 °C by a tungsten radiant heater underneath. The Nichrome 
mixture will vaporize almost as soon as it touches the tantalum strip 
and will deposit on the glass substrate. 

Separate filament evaporation requires two tungsten - filaments, one 
wrapped with nickel, the other with chromium. Both are resistance-heated 
until the two metals vaporize. This method would not give as uniform a 
film as the other because of the difference in boiling points, etc. , of the 
two metals. 

The substrate will be a very fine strand of SiO a abouj: 10 2 cm 
in diameter. As has been done in nearly all previous cases, the glass will 
be heated by a radiant heater for better adhesion and uniformity. While 
deposition is taking place, a vacuum tube voltmeter will be connected across 
the substrate to measure its resistance and to act as a control on the 
metallic evaporation. 


Ill, Possible Future Activity 

Our attention in the future will be focused on the actual con- 
struction of the resistive proportional tube. It should differ little from 
the construction of a normal proportional tube except that the strand may 
be fragile and require special consideration as far as installation. 

The first step will be the fabrication and testing of the resistive 
wire outside the tube. 
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After the assembly of the tube with several types of wires 
and variations in constructional details, the differences will be eval- 
uated and the most stable and best resistive values will be tested 
further. When the resistive proportional tube is developed to a work- 
able state, we will use it to design analytical systems for recording 
magnitude and position readings. During this stage, the tube will also 
probably evolve to a much higher degree. 


When these systems are in working order, several resistive 
proportional tubes will be linked into an array to determine whether 
they would be suitable for use in arrays. Following this, the array 
will be connected to telemetry equipment for broadcast of information. 

The next step will be to miniaturize the equipment in size 
and weight and to increase the strength and stability of the resistive 
proportional tube to withstand the rigors of space flight. Finally, 
the device will be flown in an unmanned satellite for cosmic x-ray 
observation. 
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E. Miscellaneous 


II. E. 1. 


A UNIVERSAL MAGNITUDE SYSTEM FOR ASTRONOMICAL OBJECTS 

Truman P. Kohman 

The great variety of magnitude scales in use in astronomy is confusing, 
partly because of their arbitrariness, and particularly because of the lack 
of a logical relationship between the scales for different quantum energies. 
A universal system of magnitudes based on power flux is proposed. 

This proposal is based on the concept of the expression of spectral 
(monochromatic) power flux as the logarithmic derivative with respect to any 
spectral parameter: 

I* = 


dl 


ai 


dl 


dl 

a In A 


d In v 


d In G 


d In e 


The logarithmic derivative is related to the simple derivatives by 

I* * A I' = v I* = a I* = g I* 

A v ■ 0 e 

and is proposed as a natural universal system for power fluxes. 

The universal integral and spectral magnitudes are respectively defined 

and 


by: 


- 5 , I 

m = - g log r 


5 i I* 
m* = - - log =- 


It is proposed that the scale be normalized by adopting as the fluxes for 
zero magnitude : 


I = 10 -7 * 6 
o 


W m" 2 = 2.512 X 10 " 8 ¥ m“ 2 


I* « 10 
o 


- 7. 6 


W m~ 2 un -1 = 2.512 X,10~ 8 W m~ 2 un -1 


Altemate concise formulations of the definitions of integral and spectral 
magnitudes on the universal scale become: 

5 - - 1 log Tv* - 19 ^ - I B-jS ' m -i - 19 

The system is readily extended to include partial-integral magnitudes, 
luminance magnitudes, and absolute magnitudes. 
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This system realizes' a definite physical relationship between magnitudes 
and power fluxes, facilitates meaningful intercomparisons of power fluxes or 
magnitudes in all parts of the electromagnetic spectrum from radio to gamma . 
radiation, and simplifies calculations and data presentations. 

A number of examples from various branches of astronomy illustrate the 
applicability and advantages of the system. 

This paper was presented at the meeting of the American Astronomical 
Society in Boulder, Colorado, on 1970 June 10, and draft copies (NY0-644-86) 
were given or have been sent to a number of astronomers with requests for 
comments. After consideration of comments already received and to be received, 
a final version will be submitted to the Astronomical Journal for publication. 



